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ABSTRACT 


This  document  teaches  engineers  how  to  write  computer  programs  independent  of 
any  specific  programming  language. 

It  shows  the  four  principles  of  programming  and  how  to  use  them  to  build  com¬ 
puter  programs  in  English.  It  shows  how  to  translate  the  English  programs  into 
an  appropriate  engineering  computer  programming  language:  ratfor,  Pascal,  C, 
FORTRAN77 ,  FORTRAN66,  or  BASIC. 

It  also  provides  a  programming  checklist  and  a  debugging  guide. 

ADMINISTRATIVE  INFORMATION 

This  document  was  prepared  using  Structures  Department  overhead  funds. 


INTRODUCTION 

I  presume  you  to  be  an  engineer  who  wishes  to  write  programs  to  solve  engineer¬ 
ing  problems.  You  may  already  have  had  some  experience  with  computers,  but  you 
don't  necessarily  wish  to  become  a  'computer  expert.'  Unfortunately,  use  of 
modern  computers  requires  some  dexterity  in  at  least  four  areas: 

-  you  should  have  some  acquaintance  with  your  computer's  'operating  system 
(the  program  that  lets  other  programs  run); 

-  you  must  be  able  to  wield  a  text  editor  (the  program  that  lets  you  type 
programs  and  data  into  your  computer); 

-  you  have  to  be  expert  enough  in  your  technical  field  to  know  that  you  need 
to  do  some  computing; 

-  you  have  to  express  the  solution  to  your  technical  problem  in  some  program¬ 
ming  language. 

Operating  systems  and  text  editors  have  been  addressed  in  other  documents  (see 
[Roth83l*  and  [Roth82l),  and  I  can't  begin  to  teach  you  your  field,  but  perhaps 
I  can  help  you  with  the  programming  process  itself.  ^ 

Programs  can  be  written  by  applying  a  few  simple  principles.  This  text  teaches 
these  simple  principles  in  a  language  you  already  know,  i.e.,  English.  I  will 
further  presume,  however,  that  you  are  trying  to  become  familiar  with  the  syntax 
of  one  of  the  programming  languages  illustrated  here. 

Over  85$  of  the  programs  in  the  world  are  less  than  200  lines  long,  so  you  will 
probably  spend  most  of  your  programming  time  writing  programs  of  this  size. 
Development  of  larger  programs  is  a  different  story  because  you  need  a  well- 

*  References  are  indicated  by  the  first  four  letters  of  the  author  s  name 
followed  by  the  year  of  publication,  in  brackets.  A  complete  listing  of 
references  is  given  on  page  111. 
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developed  strategy  to  manage  the  effort,  in  addition  to  knowing  the  principles 
of  programming.  In  this  slim  text,  however,  we  will  concentrate  primarily  on 
the  basic  programming  principles  (which  are  applicable  to  programs  of  any  size) 
and  postpone  discussion  of  the  management  of  large  program  development  to  a 
later  volume.  (Certainly,  you  must  be  able  to  write  a  small  program  before  you 
can  write  a  big  one.) 

The  text  is  organized  as  follows:  The  Four  Principles  of  program  construction 
are  demonstrated  in  Chapter  1 .  Chapter  2  discusses  the  passive  parts  of  com¬ 
puter  programs:  data.  Chapter  3  describes  the  Three  Parts  of  every  program. 
Chapter  4,  which  is  fairly  long,  describes  The  Constructs  (from  Control 
structures)  with  exercises  to  test  your  learning  as  you  go.  Chapter  5  describes 
the  process  of  getting  your  program  to  run,  and  Chapter  6  gives  a  checklist  to 
help  ensure  that  your  program  is  in  the  best  possible  condition  before  (and 
after)  you  run  it.  Chapter  7  is  a  short  lesson  in  correcting  your  program  when 
it  doesn't  run  correctly  (yes,  even  with  all  of  the  'good  stuff'  in  the  early 
chapters,  things  still  go  wrong).  Chapter  8  contains  capsule  reviews  of  several 
programming  languages  and  reveals  my  not-very-well-hidden  bias  for  certain 
languages.  Chapters  9  through  14  are  each  devoted  to  a  specific  programming 
language;  there  is  not  enough  detail  presented  here  to  make  full  use  of  any 
language,  but  enough  to  get  you  going  successfully.  Following  the  list  of 
references,  I  have  included  a  small  glossary  of  those  terms  having  special  jar¬ 
gon  value  to  programmers.*  I  therefore  recommend  that  you  read  Chapters  1 
through  6,  skim  7  and  8,  and  select  the  programming  language  chapter(s)  to  read 
per  your  own  taste. 

I  have  a  few  other  objectives  for  writing  this  text.  First:  algorithms  are 
often  published  in  a  language  bearing  great  resemblance  to  Pascal  and  Algol.  It 
is  desirable  to  be  able  to  read  these  algorithms,  and  perhaps  convert  them  into 
other  languages.  By  presenting  the  principles  of  programming  here  in  several 
languages,  I  hope  to  make  it  easier  for  you  to  read  and  understand  these  algo¬ 
rithms. 

Second:  (a  confession)  I  am  comfortable  using  several  kinds  of  large  mainframe 
computers  because  I  am  familiar  with  the  way  they  behave.  I  know  how  to  use 
several  text  editors  and  many  other  software  tools,  having  invented  several  of 
my  own. 

I  noticed  that  I  was  extremely  reluctant  to  program  the  Texas  Instruments  TI59 
programmable  pocket  calculator.  I  examined  my  motives  for  putting  off  working 
with  the  machine,  and  arrived  at  the  following  list  of  technical  and  emotional 
issues: 


The  user  manual  is  thick, 
crummy  program?') 


('You  mean  I  have  to  read  all  this 


write 


-  There  aren't  many  experts  around  for  me  to  talk  to  when  I  run  into  trouble, 
or,  they  would  be  busy  with  their  own  work  (which  I  would  be  disturbing);  I 
was  'on  my  own. ' 

*  If  you  don't  understand  something  on  first  reading,  skip  over  it.  Later 
text  should  make  it  clear. 
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-  I  didn't  know  what  response  the  machines  would  make  to  my  commands;  conse¬ 
quent  ’.y  I  didn't  know  what  questions  to  ask  when  something  'went  wrong.' 

-  My  program  wouldn't  work  because  I  overlooked  something  very  simple;  I 
would  therefore  appear  to  be  stupid,  and  be  laughed  at. 

-  Work  with  the  machines  is  tedious.  There  are  no  software  tools  to  help 
program  these  devices. 

-  I  don't  have  time  to  learn  how  to  use  the  machine.  And  running  speeds  are 
very  slow. 

-  I  don't  believe  programmable  pocket  calculators  are  effective,  so  I  did  not 
want  to  become  an  expert  in  programming  them. 

You  may  have  similar  apprehensions  about  using  mainframe  computers.  I  believe 
my  excuses  (the  ones  with  emotional  content)  stem  from  the  'fear  of  the 
unknown.'  So,  in  addition  to  the  reasons  mentioned  above,  I  wrote  this  to  remove 
as  many  unknowns  as  possible,  and  thereby  reduce  the  fear  level. 

Finally:  remember  that  if  you  can  afford  to  write  a  couple  of  drafts  of  a  report 
or  letter,  you  should  allow  yourself  the  'luxury'  of  doing  two  drafts  of  a  pro¬ 
gram.  Be  prepared  to  throw  the  first  draft  away;  you  will  anyway.  Relax  and 
enjoy  yourself. 


S.-' 
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1.  THE  FOUR  PROGRAMING  PRINCIPLES 

Before  you  write  any  program,  you  should  ask  yourself  three  questions: 

-  What  is  this  program  supposed  to  do? 

-  Has  someone  else  already  written  one  to  do  this? 

-  Can  several  existing  programs  be  hooked  together  to  produce  the  same 
effect? 

You  should  be  able  to  answer  the  first  question  in  one  sentence.  If  you  can  do 
what  you  want  to  do  faster  with  a  pencil  and  paper,  you  probably  don't  need  a' 
program. 

If  the  answer  to  the  second  or  third  question  is  'yes,'  you  can  often  save  your¬ 
self  a  lot  of  work  by  using  an  existing  program  instead  of  writing  a  new  one 
(you  can  also  skip  reading  this  text).  Assuming  that  a  program  doesn't  exist  to 
do  what  you  want,  or  an  existing  one  doesn't  do  it  the  way  you  want  it  done,  and 
there's  no  one  around  who  can  write  it  for  you,  you  may  decide  to  compose  your 
own.  Now  you  need  to  know  something  about  programming.* 

But  you  don't  need  to  know  a  lot!  In  fact,  you  need  only  the  following  Four 
Principles: 

-  Programs  'operate  on'  data  of  specific  types  (Chapter  2). 

-  Programs  have  at  most  three  Parts  (Chapter  3)» 

-  Programs  are  built  from  at  most  six  Constructs  (Chapter  4). 

-  Half  of  the  constructs  may  be  altered  in  at  most  two  ways  (Chapter  4). 

When  you  apply  these  principles,  you  can  compose  your  program  in  English  (or  the 
reasoning  language  in  which  you  do  your  other  thinking)  and  then  translate  the 
'English'  program  into  the  programming  language  of  choice.  The  closer  the  tar¬ 
get  language  is  to  the  way  you  think,  the  easier  the  translation  will  be. 
Often,  the  choice  of  programming  language  is  predetermined,  because  you  are 
modifying  a  program  which  already  exists,  or  you  are  writing  the  program  for 
someone  else. 

My  own  practice  is  to  think  with  a  pencil  and  paper,**  making  notes  in  a  mixture 
of  English,  rat for,  and  Pascal,  which  I  then  translate  completely  into  the  tar¬ 
get  programming  language.  Although  some  panache  is  demonstrated  by  developing 
programs  on  scraps  of  paper,  envelope  rears,  and  irregular  corners  of  slates,  it 
is  advisable  to  keep  the  notes  you  make  when  you  develop  your  program.  These 

• "Notice  that  the  initial  effort  in  writing  a  program  is  to  avoid  writing  a 
program! 

*•  Rather  than  erase  errors,  I  find  it  preferable  to  cross  them  out  and  rewrite 
correctly.  The  erossed-out  stuff  provides  a  trace  on  the  path  of  my 
thinking,  which  is  sometimes  helpful  to  have  when  debugging  time  comes. 
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will  help  you  when  your  boss  insists  that  you  document  your  work  in  a  more 
appropriate  form,  such  as  a  report.  It  is  easier,  of  course,  to  write  the 
users'  manual  for  your  program  before  you  write  the  program  itself. 

'Efficient'  programs  use  little  computer  time  and  little  computer  memory.  But: 
computer  time  and  memory  are  in  most  cases  just  too  cheap  to  worry  about!  Your 
time  is  infinitely  more  valuable  than  silicon  time.  Thus,  efficiency  is  that 
which  makes  your  life  'easier.'  Paraphrasing  Boudreau  [Boud71],  "A  big,  slow, 
correct  program  is  better  than  a  slick  wrong  one."  So  get  it  right  before  you 
make  it  'fast.'  As  you  gain  experience  using  the  four  Principles,  you  will 
notice  that  you  are  producing  programs  which  are  not  only  correct,  but  readable, 
understandable,  and  quite  fast,  as  well. 

The  action  "translating  completely  into  the  target  programming  language"  is 
called  "coding."  Because  translating  incomplete  thoughts  into  a  programming 
language  is  quite  difficult,  please  regard  the  following  maxim: 

THE  SOONER  YOU  START  CODING, 

THE  LONGER  IT  WILL  TAKE  TO  FINISH  YOUR  PROGRAM. 

When  the  translation  from  English  to  programming  language  is  complete,  the  pro¬ 
gram  can  be  typed  into  a  computer  (see  [Roth82]  for  help  in  this  area),  and  com¬ 
piled  and  executed  using  the  procedures  of  [Roth83].  Each  of  the  language 
chapters  in  this  report  also  shows  how  to  do  it. 


2.  DATA 


2.1  Data,  Constants,  and  Variables 

Data  are  the  objects  on  which  computer  programs  operate. 

Constants  are  items  of  data  which  do  not  change. 

Variables  are  things  which  may  assume  the  values  of  different  constants. 

2.2  Types 

We  can  assign  values  to,  and  examine  the  values  of,  variables  according  to  the- 

type  of  the  data. 

The  data  types  useful  for  the  majority  of  straightforward  engineering  programs 

are: 

integer  Numbers  which  have  no  fractional  part.  Most  often  used  to  count 
things.  For  example,  the  number  of  keys  on  a  piano  is  an  integer  con¬ 
stant  (88).  A  variable  used  to  store  this  data  type  is  called  an 
' integer  variable . ' 


real  Numbers  which  may  have  a  fractional  part.  Useful  for  representation 

of  physical  quantities  like  mass,  density,  thickness,  location  in 
space,  etc.  The  ratio  of  the  circumference  to  the  diameter  of  a  cir¬ 
cle  is  a  real  constant:  3.14159....  A  variable  used  to  store  this 


data  type  is  called  a  'real  variable.' 

Boolean  (or  logical).  An  item  of  Boolean  data  may  have  only  one  of  two  possi¬ 
ble  values:  'true'  or  'false.'  This  type  allows  control  of  programs 
through  logic.  It  is  'true'  that  2.0  is  greater  than  1.0.  A  variable 
used  to  store  a  Boolean  constant  is  called  a  'Boolean  variable.' 

character  Used  to  read  and  write  the  character  set  available  on  the  machine. 

These  are  typically  the  lower  case  letters,  the  upper  case  letters, 
the  digits,  the  space,  the  tab,  special  symbols  such  as  !,  @,  #,  etc. 
A  variable  used  to  store  this  data  type  is  called  a  'character  vari¬ 
able.  ' 

2.3  Expressions 

Each  of  the  types  is  defined  for  an  expression  in  a  language.  A  valid  expres¬ 
sion  is  a  combination  of  constants,  variables,  and  symbols  which  conforms  to  the 
syntax  of  the  language.  You  can  assign  the  value  of  an  expression  to  a  variable 
of  the  appropriate  type  via  the  assignment  statement  of  the  language  (more  on 
this  later). 

2.4  Data  Organizations 


The  data  organizations  most  useful  for  straightforward  engineering  programs  are 
the  simple  variable,  the  array,  and  the  file. 


DATA 


Data  Organizations 


2.4.1  The  Simple  Variable  A  simple  variable  is  a  thing  to  which  you  can  assign 
values  according  to  its  type,  compare  its  value  with  variables  or  expressions  of 
the  same  type,  and  write  the  value  of  the  variable  (so  you  can  see  the  results 
of  assignments  and  comparisons). 

All  computing  with  variables  is  done  in  accordance  with  the  variable's  type. 
For  example,  we  can  assign  "position  on  the  x-axis"  to  a  variable  we  name  X  of 
type  real.  Or  we  can  assign  the  truth  of  the  expression  "X  is  greater  than  6.1" 
to  a  variable  B  of  type  Boolean. 

2.4.2  The  Array  An  array  is  an  ordered,  determinate  number  of  variables.  You 
should  think  of  an  array  as  a  single  entity,  even  though  it  may  have  many  parts. 
For  example,  the  direction  cosines  of  a  vector  form  a  single  array  of  three  real 
data  items.  We  might  name  this  array  DIRCOS. 

The  generalization  from  simple  variable  to  array  suggests  'an  array  of  arrays,' 
'an  array  of  arrays  of  arrays,'  etc.  For  example,  we  can  consider  the  location 
of  a  point  in  space  to  be  given  by  the  three  real  coordinates  of  the  point:  X, 
Y,  and  Z  (the  names  of  three  real  3imple  variables).  Or,  we  can  consider  the 
location  to  be  given  by  an  array  which  contains  three  real  variables.  Suppose 
we  name  this  array  'X.'  Then  X(1)  could  be  the  'x'  coordinate,  X(2)  the  'y' 
coordinate,  and  X(3)  the  'z'  coordinate. 

Further,  if  we  have  10  points  in  space,  we  can  conceive  of  them  as  an  array  of 
arrays.  Suppose  we  name  the  array  COORD.  Then  we  need  10  'places'  in  the  array 
(one  for  each  point).  In  each  of  the  'places'  we  need  an  arrav  of  three  real 
variables  (one  for  each  axis).  So, 

COORD ( 1,1)  might  be  the  'x'  coordinate  of  the  first  point, 

COORD (4, 2)  could  be  the  'y'  coordinate  of  the  fourth  point, 

COORD (7, 3)  could  be  the  'z'  coordinate  of  the  seventh  point. 

In  the  COORD  array,  therefore,  the  first  'subscript'  gives  the  number  of  the 
point  and  the  second  gives  the  axis  of  measurement.  The  important  thing  to  note 
is  that  the  array  concept  allows  us  to  deal  with  a  lot  of  numbers  with  a  single 
name . 

Note  that  although  each  item  of  an  array  can  be  of  any  type  (the  example  above 
used  an  array  of  type  real) ,  the  subscripts  are  always  of  type  integer.  Ry 
using  integers  in  control  structures,  we  can  process  arrays  very  quickly  and 
with  great  economy  of  language. 

2.4.3  The  File  A  file  is  an  ordered  collection  of  variables  which  may  be  of 
unknown  length.  To  indicate  a  file,  I  will  use  a  name  surrounded  by  brokets: 
<input>  is  to  be  read  'the  file  whose  name  is  input,'  or  simply  'the  input 
file.'  Examples: 

-  <input>  to  a  finite  element  program  is  an  unknown  number  of  'lines,'  each 
line  being  an  array  of  80  variables  of  type  character. 
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-  This  report  originated  as  several  files  on  a  UNIX*  system  (at  the  time  I 
first  wrote  this  sentence,  the  files  were  2306  lines  long,  and  had  9020 
'words'  made  up  of  51560  characters). 

Files  are  used  for  input  and  output  to  computer  programs.  Although  the  'memory' 
of  a  computer  is  indeed  finite,  the  files  available  to  a  program  are  conceptu¬ 
ally  infinite.  This  allows  us  a  handy  place  to  temporarily  save  large  numbers 
of  simple  variables  and  arrays  which  would  exceed  the  capacity  of  the  computer's 
memory. 

Because  of  the  possibility  that  files  can  be  very  long,  programming  languages 
usually  provide  access  to  files  through  a  "window; "  that  is,  you  only  get  to  see 
a  little  of  it  at  a  time.  The  shape  of  this  window  is  greatly  dependent  on  the 
operating  system  and  the  language  being  used. 

2.4.4  Other  types  and  organizations  Although  the  types  and  organizations  men¬ 
tioned  above  are  satisfactory  for  the  vast  majority  of  the  programs  which  you 
will  write,  the  list  is  incomplete.  FORTRAN,  for  example,  provides  the  addi¬ 
tional  types 

complex  A  variable  which  has  a  'real'  and  'imaginary'  part;  both  of  the  parts 
are  represented  internally  as  'reals'. 

double  A  variable  which  has  twice  the  precision  of  a  'real.' 

Other  languages  provide  other  types  and  organizations  which  include  the  SET,  the 
RECORD,  the  TREE,  etc.  To  use  these  in  a  language  which  does  not  explicitly 
provide  them  requires  ingenuity  in  mapping  the  structures.  You  are  limited  only 
by  the  limits  of  your  imagination,  and  the  capability  to  express  your  idea  in 
terms  of  a  known  language  running  on  an  available  machine. 

The  best  treatment  of  data  types  and  organizations  is  given  by  [Wirt76],  which 
is  based  on  the  Pascal  language.  In  FORTRAN,  [Day72]  and  [Berz71]  are  quite 
good. 


*  UNIX  is  a  trademark  of  Bell  Laboratories. 
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2.5  Variable  Names 

"Don't  stand  there  chattering  to  yourself  like  that," 

Humpty  Durapty  said,  looking  at  her  for  the  first  time,  "but 
tell  me  your  name  and  your  business." 

"My  name  is  Alice,  but — " 

"It's  a  stupid  name  enough!"  Humpty  Dumpty  interrupted 
impatiently.  "What  does  it  mean?" 

"Must  a  name  mean  something?"  Alice  asked  doubtfully. 

"Of  course  it  must,"  Humpty  Dumpty  said  with  a  short 
laugh:  "my  name  means  the  shape  I  am  -  and  a  good  handsome 
shape  it  is,  too.  With  a  name  like  yours,  you  might  be  any 
shape,  almost."  [Carr96] 

Taking  the  Humpty  Dumpty  hint,  we  will  wish  to  name  variables  in  such  a  way  that 
their  names  give  us  a  clue  to  their  essence.  Because  English  is  our  language, 
we  have  no  limitations  on  the  length  or  style  of  a  variable  name.  What  is 
interesting  about  English  (or  mathematical)  variable  names  is  that  they  are 
assigned  in  a  context  which  preserves  their  meaning.  Thus,  when  we  say  some¬ 
thing  like,  'Let  P  be  a  point  in  space,'  then  the  name  'P'  somehow  is  seen  to  be 
tied  to  the  word  'point.'  Of  course,  a  second  point  would  need  another  name  to 
make  it  distinct  from  P.  Mathematicians  will  usually  try  to  pick  a  name  that  is 
somehow  close  to  the  first  name,  and  yet  obviously  different;  perhaps  'Q'  (it  is 
'near'  P  alDhabetically) .  We  note  that  P  and  Q  are  short  names  because  we  usu¬ 
ally  prefer  not  to  write  very  much.  At  the  same  time,  names  (i.e.,  words)  which 
are  longer  than  eight  characters  are  usually  at  least  three  syllables,  which 
makes  them  hard  to  read  and  harder  to  type. 

On  the  other  hand:  a  textbook  which  defines  'P'  on  page  2,  say,  and  never 
reminds  us  that  P  is  a  point,  is  going  to  be  a  difficult  text  to  refer  to. 
That's  why  mathematicians  are  always  redefining  things  for  their  readers;  con¬ 
text  is  usually  limited  to  about  a  page  or  two.  You  should  cultivate  the  same 
habit  in  your  programs,  because  programs  are  not  always  read  from  beginning  to 
end;  one  enters  them  at  any  place  suitable  to  the  need  at  hand.  It  is  therefore 
important  that  the  program  reader  be  able  to  get  a  fair  understanding  of  what  is 
going  on  from  the  immediate  context. 

Thus: 

VARIABLE  NAMES  SHOULD 
BE  OF  'REASONABLE'  LENGTH. 

'Reasonable'  means  from  one  to  eight  characters  in  length  (maximum  length  of 
twelve  characters),  or  redefined  with  commentary  at  fairly  frequent  intervals  in 
the  program  text.  In  the  'Let  P  be  a  point  in  space'  case,  my  practice  would  be 
to  name  the  variable  P  in  a  program  of  less  than  50  lines,  and  name  it  POINT  in 
any  program  of  50  lines  or  larger. 

The  creative  act  of  naming  variables  is  difficult!  A  thesaurus  can  be  quite 
helpful. 
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Exercise:  Obtain  someone's  program  (preferably  one  of  your  own  that  la  at  least 
6  months  old),  and  take  a  look  at  the  variable  names.  Can  you  tell  what  they 
are?  Why  or  why  not? 


3 .  THE  THREE  PARTS  OF  EVERY  PROGRAM 


Computer  programs  have  three  parts: 

-  The  Initialization. 

-  The  Computation. 

-  The  Clean-up. 

3.1  The  Initialization 

This  program  part  is  that  in  which  we  get  ready  to  compute.  It  includes 
the  program  header  line, 
the  declaration  of  variable  data  types, 
the  opening  and  positioning  of  files, 
the  setting  of  default  values,  and 
the  setting  of  initial  conditions. 

On  some  computers,  the  operating  system  itself  provides  much  of  the  program  ini¬ 
tialization.  The  pieces  most  often  under  your  control  are  the  declaration  of 
variable  type,  and  the  initialization  of  variables  to  some  value .  The  declara¬ 
tions  are  much  like  the  statements  preceding  a  mathematical  proof  ("Let  f  be  a 
smooth  continuous  function  in  the  domain  D  with  f '  the  first,  and  f"  the 
second,  derivative  of  f  with  respect  to  distance  in  D”). 

Initialization  of  a  variable  to  a  specific  value  uses  the  'assignment  state¬ 
ment,'  which  'assigns'  a  value  to  a  variable.  Because  the  form  of  the  assign¬ 
ment  statement  varies  from  language  to  language,  each  language  chapter  includes 
a  section  defining  this  statement. 

The  initialization  is  usually  the  last  part  of  a  program  to  be  completed.*  This 
i3  because  you  don't  know  what  the  initial  values  should  be  when  you  begin  com¬ 
position  of  the  program.  Hence,  initialization  can  also  be  the  hardest  part  of 
a  program  to  write. 


*  This  is  comparable  to  the  way  in  which  a  report  or  book  is  written:  the  last 
parts  to  be  written  are  the  introduction  and  the  abstract. 
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3.2  The  Computation 


THE  THREE  PARTS  OF  EVERY  PROGRAM 


This  is  the  program  part  which  does  your  work.  It  is  the  part  where  you  will 
use  The  Constructs. 

3.3  The  Clean-up 

In  this  part  of  the  program,  we  save  those  things  which  we  need  to  keep,  throw 
away  the  garbage,  and  exit  from  our  computations  gracefully  (if  possible). 

On  many  computers,  the  entire  clean-up  is  performed  by  the  operating  system. 
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4.  THE  CONSTRUCTS 


In  this  fairly  long  chapter,  I  present  The  Constructs  both  in  English  and  graph¬ 
ically. 

4.1  Boolean  Logic 

The  facility  for  logic  built  into  programming  languages  is  a  mechanical  imple¬ 
mentation  of  Boolean  logic.  Boolean  expressions  are  those  which  allow  us  to 
test  truth  and  falsity  according  to  the  rules  of  logic  formalized  by  Boole.  A 
Boolean  expression  may  have  only  one  of  two  values:  TRUE  or  FALSE. 

4.2  Form  of  Boolean  Expressions 


In  the  following,  X  and  Y  are  variables  or  expressions  of  the  same  type 
each  'Boolean  expression'  will  be  either  TRUE  or  FALSE. 

X  is  greater  than  Y. 

X  is  greater  than  or  equal  to  Y. 

X  is  equal  to  Y. 

X  is  not  equal  to  Y. 

X  is  less  than  Y. 


X  is  less  than  or  equal  to  Y. 


,1 
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The  actual  activity  of  a  program  is  performed  by  statements  which  conform  to  the 
language  syntax.  A  statement  in  a  language  will  be  represented  here  in  the  form 
'S',  'SI',  ' S2 ' ,  etc. 

One  of  the  most  important  statement  types  is  the  assignment  statement,  which  has 
the  form 

X  ao  E 

in  which  'X'  is  a  variable,  'ao'  is  the  assignment  operator  (usually  a  symbol 
composed  of  1  or  more  characters),  and  'E'  is  an  expression  of  the  same  type  as 
'X'.  One  reads  an  assignment  statement:  "the  current  value  of  'X'  is  replaced 
with  the  value  'E'"  or,  more  simply,  "'X'  becomes  'E'."  This  is  obviously  dif¬ 
ferent  from  the  statement  "'X'  equals  'E'.w  The  latter  is  a  'Boolean  expression' 
which  may  assume  only  the  value  TRUE  or  FALSE;  the  former  is  called  the  'assign¬ 
ment  statement'  because  the  value  of  an  expression  is  assigned  to  a  variable.* 
Note:  execution  of  the  assignment  statement  destroys  any  previous  value  of  X! 

We  emphasize  that  the  assignment  statement  is  not  the  same  as  the  Boolean 
expression  'X  equals  E';  the  assignment  statement  is  an  action,  whereas  the 
Boolean  expression  is  an  assertion. 

For  example,  let's  assume  that  the  integer  variable  X  has  the  value  1.  Then 
X  equals  X  +  1 
is  false,  whereas 
X  becomes  X  +  1 
assigns  the  value  2  to  X. 

A  statement  will  be  diagramed 


where  the  arrow  into  the  box  Indicates  that  control  is  passing  to  the  statement 
'S.'  When  the  statement  'S'  has  completed  its  action,  control  passes  to  the 
'next'  statement  to  be  executed;  this  is  represented  by  the  arrow  out  of  the 
box. 

*  The  expression  may  be  as  simple  as  another  variable,  in  which  case  we  say 
that  one  variable  is  copied  into  another. 


-14- 


THE  CONSTRUCTS 


The' Assignment  Statement 


In  general ,  the  notation  'S'  may  be  read  as  'n  statements  in  sequence'  where  "n" 
is  greater  than  or  equal  to  0.  Thus, 


may  actually  represent  something  like 


Most  contemporary  computers  are  sequential  machines,  which  means  that  statements 
are  executed  in  the  order  in  which  they  are  accessed.  Thus,  if  we  write  the 
statements  SI,  followed  by  S2,  followed  by  S3,  the  statements  will  be  executed 
in  the  order  SI,  S2,  S3  unless  we  somehow  change  this  order.  This  'changing'  is 
done  with  The  Constructs.  Please  note  that  each  of  The  Constructs  is  also  a 
' statement j'  and  may  be  hooked  together  Just  as  anv  of  the  'other'  statements  of 
a  language  are. 
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The  Conditional  Construct 

The  Conditional  construct  permits  computations  to  be  performed  only  if  certain 
conditions  are  met.  Given  that  B  is  a  Boolean  expression,  the  Conditional  con¬ 
struct  has  the  form 

If  B, 

then  S. 

For  example: 

If  "it  is  raining" 

then  "carry  an  umbrella." 

In  this  example,  the  Boolean  expression  B  is:  "it  is  raining."  This  expression 
may  be  TRUE  or  FALSE.  The  action  to  be  performed,  S,  is  "carry  an  umbrella." 

The  Conditional  construct  may  be  diagramed 


where  the  'if'  part  of  the  statement  is  represented  by  the  B  inside  the  diamond. 
If  'B'  is  true,  then  the  true  path  is  taken,  and  'S'  is  done.  If  'B'  is  false, 
then  the  false  path  is  taken. 
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The  Alternative  Construct 


When  there  are  alternative  actions  which  may.  be  performed,  we  can  use  the  Alter¬ 
native  construct.  Given  that  B  is  a  Boolean  expression,  then  the  Alternative 
construct  has  the  form 

If  B 

then  SI ; 
else, 

S2. 

For  example: 

If  "it  is  raining" 

then  "carry  an  umbrella;" 
else, 

"lower  convertible  top." 

Note  that  these  actions  are  mutually  exclusive.  That  is,  "if  it  is  raining" 
then  the  only  statement  which  is  executed  is  "carry  an  umbrella."  "If  it  is  not 
raining,”  then  the  only  statement  which  is  executed  is  "lower  convertible  top." 

The  Alternative  construct  can  be  diagramed 
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The  Multiple  Choice  Construct  THE  CONSTRUCTS 

4.6  The  Multiple  Choice  Construct 

Out  of  a  number  of  possibilities,  we  may  wish  to  select  the  appropriate  one  and 
do  the  statement(s)  associated  therewith.  For  example: 

In  case  of 

rain:  "carry  an  umbrella" 
snow:  "wear  parka  and  hat" 
hail:  "stay  home" 
tornado:  "goto  cellar" 
sunshine:  "play  tennis" 

where  only  one  of  the  choices  is  elected.  Thus,  "goto  cellar"  and  "play  tennis"' 
are  mutually  exclusive  options. 

Another  way  of  expressing  this  construct  is  with  a  succession  of  Alternative 
constructs: 

if  "it  is  raining"  then 
"carry  an  umbrella" 
else  if  "it  is  snowing"  then 
"wear  parka  and  hat" 
else  if  "it  is  hailing"  then 
"stay  home" 

else  if  "it  is  tornadoing"  then 
"goto  cellar" 

else  if  "it  is  sunshining"  then 
"play  tennis" 
else 

"turn  on  radio" 


where  we  have  included  the  last  'else'  to  take  care  of  the  time  when  none  of  the 
weather  conditions  we  know  about  occurs. 


Iteration 
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4.7  Iteration 


Iteration  is  one  of  the  great  strengths  of  computers;  the  machines  will  do  the 
same  actions  over  and  over  with  no  complaint. 

There  are  two  types  of  iteration:  indeterminate,  and  determinate.  Indeterminate 
iteration  can  be  profitably  applied  to  reading  numbers  from  files  when  the 
length  of  the  file  is  unknown  prior  to  reading.  Determinate  iteration  may  be 
used  when  certain  portions  of  an  array  whose  size  is  completely  known  must  be 
set  to  some  value. 

Indeterminate  iteration  may  be  classified  into  two  subtypes: 

-  Test,  then  (perhaps)  act.  This  iteration  construct  is  called  the  While. 

-  Act,  then  test.  This  construct  is  called  the  Repeat-Until. 
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4.7.1  The  While  Construct  The  indeterminate  iteration  is  applicable  in  those 
cases  when  the  number  of  times  the  iteration  will  be  done  is  not  known  a  priori. 
Given  that  B  is  a  Boolean  expression,  then  the  While  construct  has  the  form 

while  B, 

S. 

Note  that  the  statement  S  may  not  be  executed  at  all!  For  example: 

while  "not  quitting  time" 

"do  some  more  work" 

If  it  is  quitting  time,  then  "do  some  more  work"  is  not  done.  Another  example: 

While  "you  have  not  reached  the  end  of  a  file" 

"read  and  process  the  next  line  of  the  file" 

If  we  are  already  at  the  end  of  file,  then  "read  and  process  the  next  line  of 
the  file"  is  not  even  attempted. 

The  While  construct  can  be  diagramed 


The  'going  back'-ness  of  this  construct  is  what  gives  rise  to  the  term  'loop'  to 
describe  iteration.  I  will  henceforth  use  the  terms  'iteration  construct'  and 
'loop'  interchangeably. 
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4.7.2  The  Repeat-Until  Construct  This  iteration  places  the  test  at  the  end  of 
the  loop.  Note  the  difference  between  the  While  and  the  Repeat-Until:  the 
statements  controlled  by  the  While  may  not  be  executed,  whereas  the  statements 
under  control  of  the  Repeat-Until  are  guaranteed  to  be  executed  at  least  once. 
Given  that  B  is  a  Boolean  expression,  then  the  Repeat-Until  construct  has  the 
form 

Repeat 

51. 

52. 

Until  B. 

In  this  construct,  the  statements  SI  and  S2  are  guaranteed  to  be  executed  at 
least  once.  For  example: 

Repeat 

"brush  uppers" 

"brush  lowers" 

Until  "teeth  are  clean" 

The  Repeat-Until  can  be  diagramed 
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4.7.3  The  Determinate  Iteration  Construct  This  kind  of  iteration  is  called 
determinate  because  the  number  of  times  the  loop  will  be  executed  is  dependent 
on  an  index  to  some  countable  sequence.  The  .sequence  is  assumed  to  begin  at  J 
and  continue  to  K  in  increments  of  M  (the  sequence  may  increase  or  decrease). 

Let  B  be  a  Boolean  expression  which  compares  the  values  of  the  integer  variables 
I  and  K  (say,  for  example,  B  is  the  expression  "I  is  less  than  K").  Then  the 
determinate  iteration  construct  has  the  form 

I  becomes  J. 

While  B  { 

S. 

I  becomes  I  +  M. 

} 

where  the  braces  ' { '  and  ' } '  define  the  complete  range  of  statements  under  con¬ 
trol  of  the  While.  Note  that  there  exists  the  possibility  that  statements  "S" 
and  "I  becomes  I  +  M"  may  not  be  executed  at  all  (why?) . 

Determinate  iteration  may  be  diagramed 


where  the  statement  II  represents  the  initialization  of  the  counter,  and  12 
represents  the  resetting  of  the  counter  so  it  points  to  the  next  item  of 
interest. 

As  an  example  of  this  construct,  let's  find  and  print  the  largest  element  MAX  in 
1  an  array  of  integers  ZOT  which  has  N  elements.  Let's  say  N  =  3  and  the  contents 

of  ZOT  are 

ZOT( 1 )  =  27 
Z0T(2)  =  31 
ZOT( 3)  =  12 

Let's  assume  we  have  an  integer  variable  'I'  which  we  can  use  to  access  each  of 
the  elements  of  the  array  ZOT.  We'll  also  assume  that  the  contents  of  ZOT  are 
calculated  in  some  other  part  of  our  program,  and  that  we  don't  know  what  the 
values  of  ZOT  are.  We  do  know  that  there  are  N  values  we  have  to  look  at,  and 
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that  we'd  better  look  at  them  all.  So,  here  is  an  English  'code  fragment'  which 
will  allow  us  to  find  MAX 

MAX  becomes  Z0T(1). 

I  becomes  1 . 

While  "I  is  less  than  or  equal  to  N"  { 

If  "MAX  is  less  than  ZOT(I)"  then  "MAX  becomes  ZOT(I)" 

I  becomes  1+1 . 

} 

Print:  "The  largest  value  in  the  ZOT  array  is  "  MAX. 

In  this  example,  we've  combined  the  While  and  the  Conditional  constructs. 

Exercise:  mark  the  three  parts  of  this  program  fragment  to  show  the  initializa¬ 
tion,  computation,  and  clean-up  phases. 

Now  let's  step  through  the  calculations  by  substituting  numerical  values  for 
each  variable: 

MAX  becomes  27. 

I  becomes  1 . 

While  "1  is  less  than  or  equal  to  3"  { 

If  "27  is  less  than  27"  then  (not  true,  so  do  nothing!) 

I  becomes  2 

While  "2  is  less  than  or  equal  to  3"  { 

If  "27  is  less  than  31"  then  "MAX  becomes  31." 

I  becomes  3 

While  "3  is  less  than  or  equal  to  3"  { 

If  "31  is  less  than  12"  then  (not  true,  do  nothing) 

I  becomes  4 

While  "4  is  les3  than  or  equal  to  3"  (not  true,  end  the  While) 

The  largest  value  in  the  ZOT  array  is  31.  (printed  answer) 

and  the  loop  terminates  with  MAX  =31  and  1=4. 

Exercises: 

1.  Initialize  I  to  2  instead  of  1 .  Is  the  program  any  faster?  By  how  much? 
Is  it  easier  or  harder  to  understand?  Why? 

2.  Does  the  changed  program  work  if  N  =  1? 

3.  How  would  you  change  this  program  to  make  it  work  for  N  =  0?  What  value 
would  you  assign  to  MAX? 


4.8  Altered  Loops 


Strictly  speaking,  we  can  structure  our  code,  with  only  the  Conditional,  the 
Alternative,  the  Multiple  Choice,  the  While,  the  Repeat-Until,  and  the  Deter¬ 
minate  Iteration  constructs.  Using  only  these  constructs  sometimes  makes  for 
extremely  convoluted  programs,  and  the  meaning  of  a  program  can  actually  be  made 
clearer  by  performing  a  test  inside  a  loop  and  directing  the  path  of  the  execu¬ 
tion  according  to  the  test  results.  This  kind  of  construct  thus  alters  the 
nature  of  the  loop  to  'semi-determinate.'  Be  this  as  it  may,  please  observe  the 
following  maxims 

ALTER  LOOPS  INFREQUENTLY. 

You  will  notice  in  the  graphics  for  altered  loops  that  there  are  seven  paths 
between  parts  of  the  construct.  Although  this  is  still  within  the  realm  of 
comprehension,  it  is  complicated,  and  makes  understanding  of  a  program  more  dif¬ 
ficult. 
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The  first  unusual  circumstance  is  'stop  iterating  immediately  (i.e.,  leave  the 
loop);'  the  next  statement  to  be  executed  is  the  one  which  follows  the  termina¬ 
tor  of  the  loop.  This  construction  is  typical  in  parts  of  programs  which  check 
data  for  errors  and  reach  an  impasse. 

While  B1  { 

SI. 

If  B2 

then  "break" 

Else, 

S2. 

} 

S3. 

In  this  construct,  the  word  "break"  means  that  computation  is  to  continue  with 
statement  S3,  which  is  entirely  outside  of  the  range  of  the  While. 

The  Loop  Exit  can  be  used  to  alter  any  iteration  construct;  it  is  diagramed 
modifying  a  While. 
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Altered  Loops 


4.8.2  Loop  Redo 

The  second  unusual  circumstance  is  'immediately  continue  iteration  with  the  next 
test. ' 

While  B1  { 

31. 

If  B2 

then  "next" 

Else, 

S2. 

} 

S3. 

In  this  construct,  the  word  "next"  means  that  computation  is  to  proceed  with  the 
next  test  B1.  That  is,  "next"  says  "continue  computation  with  the  statement 
'While  B1'."  By  the  way,  use  of  the  "next"  is  extremely  rare  (I  have  used  it 
infrequently,  and  find  few  uses  in  the  references). 

The  Loop  Redo  can  be  used  to  alter  any  iteration  construct;  it  is  diagramed 
modifying  a  While. 
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Cy 

4.9  Construct  Summary 

Conditional: 

If  B  then 

o*- 

S. 

i 

Alternative: 

If  B1  then 

8 

kV 

SI. 

Else 

S2. 

1 

Multiple  Choice: 

In  case  of 

While: 


Repeat-Until: 


Determinate: 


Loop  exit: 
Loop  Redo: 


B1:  SI. 
B2:  S2. 
B3:  S3. 
B4:  S4. 


(or) 


If  B1  then 

51. 

Else  if  B2  then 

52. 

Else  if  B3  then 

53. 

Else  if  B4  then 

54. 


While  B 
S. 

Repeat 

S 

Until  B. 

I  becomes  J. 

While  B  { 

S. 

I  becomes  I  +  M. 

} 

Break. 

Next. 
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Which  Construct? 


Which  Construct? 


We  now  have  a  complete  set  of  constructs  to  apply  to  any  programming  problem. 
The  nagging  question  remains:  Which  construct  should  be  used  where?  You  will 
usually  have  little  difficulty  in  selecting  from  the  Conditional,  Alternative, 
or  Multiple  Choice  constructs.  The  problem  usually  lies  with  the  loops. 

The  following  may  help  in  selecting  the  appropriate  iteration  construct: 

-  How  the  iteration  is  supposed  to  stop  will  hint  at  the  form  of  the  con¬ 
struct.  By  the  way,  make  sure  it  will  stop! 

-  You  can  always  use  a  While  (you  might  have  noticed  that  the  two  other 
iteration  constructs  are  really  variations  of  the  While). 

-  If  you  can  count  the  data  items  to  be  processed,  then  the  Determinate 
iteration  is  usually  proper. 
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4.11  Examples 

These  examples  are  a  good  opportunity  for  some  exercise.  For  each  case,  write 
an  'English'  program,  and  compare  it  to  the -sample  solution  shown  at  the  end  of 
the  chapter.  Then  translate  your  program  into  a  programming  language  and  com¬ 
pare  your  translation  with  the  example  shown  in  the  section  devoted  to  that  tar¬ 
get  language. 

4.11.1  Temperature  Conversion  Produce  a  table  of  equivalent  Celsius  tempera¬ 
tures  for  the  Fahrenheit  temperatures  from  -40F  to  100F  in  increments  of  5F. 

Given  any  Fahrenheit  temperature  F  (i.e.,  F  is  the  REAL  temperature),  the  Cel¬ 
sius,  or  PHONY,  temperature  C  is  given  by  the  equation 

C  =  5*(F  -  32)/9 

where  the  symbol  '•'  means  'multiply,'  and  the  symbol  '/'  means  'divide.' 

We  might  proceed  as  follows.  Beginning  at  the  top  of  a  page  of  lined  paper, 

write  down  the  first  Fahrenheit  temperature  to  be  converted:  -40.  We  now  plug 

-40  into  our  recipe  for  Celsius  temperature 

5*(-40-32)/9  =  5«(-72)/9  =  -360/9  =  -40 

and  write  down  the  answer  -40C  in  a  second  column.  Since  we  haven't  completed 
our  job  yet  (we  haven't  reached  100F),  we  add  5F  to  the  current  temperature 
(-40F)  to  obtain  -35F.  With  Fahrenheit  temperature  of  -35F,  we  plug  and  chug, 

producing  -31C.  We  continue  this  process  until  the  table  is  full. 

Space  is  provided  here  for  your  program. 
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4.11.2  Nearest  Points  Given  10  points  in  3-D  space,  find  the  two  points  which 
are  nearest  to  each  other. 

The  solution  to  this  problem  is  intuitively  obvious.  We  merely  compute  the  dis¬ 
tance  from  each  point  to  its  neighbors  using  the  distance  formula 

d  =  sqrt[(x2-x1)»»2  +  (y2-y1)#*2  ♦  {zl-zy)**^ 

where  'd'  is  the  distance,  'sqrt'  represents  the  square  root,  '••2'  represents 
squaring  a  number,  'x2'  is  the  x-eoordinate  of  "one"  point  (y2  and  z2  are  the  y- 
and  z-coordinates  of  this  same  point),  'xl'  is  the  x-coordinate  of  "another" 
point  (yl  and  zl  are  the  y-  and  z-coordinates  of  this  same  point). 

Those  two  points  which  are  at  the  minimum  distance  are  the  two  closest  neigh¬ 
bors. 

Additional  exercises: 

1.  Modify  your  program  to  find  the  points  which  are  farthest  apart. 

2.  Print  all  pairs  of  points  which  are  at  the  same  minimum  (maximum)  distance. 

3.  Print  a  table  of  the  distances  from  each  point  to  its  neighbor. 
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4.11.3  Count  the  'A's  Read  <input>  and  count  the  number  of  'A's  in  the  file. 


This  is  a  simple  problem  mentally,  but  is  quite  tedious  and  prone  to  error. 
Let's  assume  that  we  can  tell  when  we  have  reached  the  end  of  a  file,  much  like 
we  can  tell  when  we've  reached  the  end  of  a  book.  Then  we  need  merely  read  the 
file  a  character  at  a  time,  and  check  to  see  if  the  character  we've  just  read  is 
an  'A.'  If  it  is,  then  we  can  increase  our  count  of  'A's  by  one.  (Hint:  What 
does  the  initial  value  of  our  count  have  to  be,  i.e.,  before  we  read  any  charac¬ 
ters?  Does  your  program  give  the  correct  count  when  there  are  no  'A's?  When 
there  are  no  characters  at  all?) 

Additional  exercise:  modify  your  program  to  read  real  numbers  from  a  file  and 
count  the  values  greater  than  some  threshold,  say,  10.0.  Assume  that  the  last 
!  line  in  the  file  will  have  the  number  -99999.0. 
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Solutions  to  the  Examples 


4.12  Solutions  to  the  Examples 
4.12.1  Temperature  Conversion 
Initialization  -  declarations: 

Let  F  and  C  be  real  variables  which  measure  temperature. 
Initialization  -  setting  initial  values: 

F  becomes  -40. 

Computation: 

While  F  is  less  than  or  equal  to  100  { 

C  becomes  5*(F-32)/9. 

Print:  F,  C. 

F  becomes  F  +  5. 

} 

Cleanup: 


( none ! ) 
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4.12.2  Nearest  Points 
Initialization  -  declarations: 

Let  OLDD  be  a  real  variable  measuring  the  smallest  distance. 

Let  NEWD  be  a  real  variable  measuring 
the  distance  between  any  two  points. 

Let  X  be  an  array  of  reals  containing  the  3  coordinates  of 

each  of  the  10  points.  Let  the  first  subscript  designate  the  point, 
and  the  second  designate  the  coordinate.  X  is  therefore  10  by  3. 

Let  I,  J,  K,  M,  and  N  be  integer  variables  used  for  counting. 

Initialization  -  setting  initial  values: 

Read  all  the  coordinates  into  X  from  <input>. 

OLDD  becomes  the  distance  between  point  1  and  point  2  (why?); 
i.e.,  OLDD  becomes  SQRT( (X(1,1)-X(2,1) )**2 

+  (X(1,2)-X(2,2) )**2 
+  (X( 1 ,3)-X(2, 3) )#*2) 


Computation: 

I  becomes  1 . 

While  "I  is  less  than  or  equal  to  9"  { 

J  becomes  1+1 . 

While  "J  is  less  than  or  equal  to  10"  { 

NEWD  becomes  SQRT( (X(I, 1 )-X( J, 1 ) )**2 
+  (X(I,2)-X( J,2))**2 
+  (X(I,3)-X(J,3))**2) 


(found  a  new  nearest?) 

If  "NEWD  is  less  than  or  equal  to  OLDD"  then  { 

M  becomes  I. 

N  becomes  J. 

OLDD  becomes  NEWD. 

} 

J  becomes  J+1 . 

}  (end  of  While  "J  is  less...) 

I  becomes  1+1 . 

}  (end  of  While  "I  is  less...) 

Cleanup: 

Print:  '  Closest  points  are:  ' 

Print:  X(M, 1 ) ,X(M,2) ,X(M,3) 

Print:  X(N, 1 ) ,X(N,2) ,X(N,3) 
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4.12.3  Count  the  'A's 

Initialization  -  declarations: 

Let  Count  be  an  integer  variable. 

Let  Ch  be  a  character  variable. 

Initialization  -  setting  initial  values: 

Count  becomes  0. 

Read  Ch. 

Computations: 

While  "not  end  of  report"  { 

If  Ch  =  'A'  then 
Count  becomes  Count  +  1 . 

Read  Ch.  (destroying  what  was  in  Ch) 

} 

Cleanup: 


Print:  Count 


5.  GETTING  YOUR  PROGRAM  TO  RUN 


Assuming  that  you've  written  your  program  correctly  in  'English'  and  correctly 
translated  it  into  a  programming  language,  how  do  you  get  it  onto  the  machine 
and  running? 


Presuming  you  are  working  with  a  CDC  machine,  the  steps  necessary  are: 

-  typing  the  program  into  the  computer  using  a  text  editor  (with  ED, 
example  [Roth82]), 

-  'compiling'  the  code  into  'machine  language,' 

-  'loading'  the  machine  language, 

-  and  'executing'  the  program. 

A  sketch  of  this  procedure  is  shown  below. 


ED  , 

i 

(coda) 


(input) 


COMPILER 

I""-— 

Ogo) 

I 

LOADER 

I"-—"'  / 

(■bi)  / 

I 

►  SYSTEM  - ►  (output) 


*  / 
*  t 
"  / 


-►  ERROR 

s* 

s  / 

/  / 


for 


I've  called  the  file  of  code  that  you  write  with  ED  '<code>.'  It  will  be  a  file 
of,  say,  FORTRAN,  rat for,  or  Pascal. 

A  COMPILER  is  a  program  which  translates  a  'higher  order  language'  (your  code) 
into  'machine  language,'  which  is  very  difficult  to  read  and  understand,  and 
which  varies  from  one  brand  of  machine  to  another.  Compilers  were  invented  to 
reduce  your  need  to  look  at  the  machine  language,  and  express  yourself  in  a  way 
which  is  portable  from  one  machine  to  another. 

The  machine  language  file  is  called  <lgo>  (from  load  and  go).  It  Is  difficult 
(impossible)  to  read  <lgo>  with  ED;  if  you  do~ enter  ED  with  a  <lgo>  file,  you 
will  destroy  it  if  you  do  a  SAVE! 
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'Loading'  is  a  process  which  allows  you  to  use  things  which  other  people  have 
written,  including  the  transcendental  functions,  the  input/output  routines, 
graphics  routines,  etc.  Only  files  of  machine  language  are  LOADed,  so  the  pro¬ 
cess  is  largely  invisible  to  a  casual  programmer.  Naturally  enough,  LOADing  is 
performed  by  a  program  called  a  LOADER,  which  reads  your  <lgo>  file  and  writes 
its  product,  the  'absolute'  program,  into  a  file  called  <abs>.  <abs>  is  the 
only  kind  of  file  which  can  execute  on  a  CDC  machine  (<abs>  is  also  impossible 
to  read  with  ED). 

Execution  is  the  performance  of  work  according  to  the  program  you  wrote.  This 
is  done  by  typing  the  command  'abs.' 

The  solid  lines  in  the  sketch  show  the  path  of  intellect  when  everything  'goes- 
right.'  Unfortunately,  everything  going  right  on  the  first  shot  is  a  rare 
occurrence.  Still,  you  want  to  spend  as  little  time  as  necessary  in  this  'loop' 
trying  to  get  your  program  to  run.  You  can  use  the  aids  in  the  next  two 
chapters  to  make  everything  'go  right'  as  soon  as  possible. 

The  dotted  arrows  show  the  path  of  intellect  when  something  'goes  wrong.'  You 
eliminate  errors  by  repairing  <eode>  and  repeating  the  compile-load-go  process 
until  no  errors  exist,  since  the  tools  used  here  are  very  reliable.  That  is, 
except  in  extremely  rare  cases,  the  errors  are  yours ,  not  the  computer  systems'. 

The  easiest  ways  to  effect  this  loop  quickly  are  through  the  use  of  the  pro¬ 
cedures  in  Chapter  of  [Roth83]. 
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6.  THE  PROGRAM  CHECKLIST 


Although  programs  can  be  designed  and  coded  from  simple  principles,  the  devil 
hides  In  the  details.  To  prevent  this  ancient  foe  from  getting  a  toehold  on 
their  code,  experienced  programmers  can  often  be  observed  practicing  seemingly 
weird  and  lengthy  incantations  before  they  run  their  programs.  The  list 
presented  in  this  chapter  extends  one  found  in  [Myer79]  and  initiates  the  novice 
into  this  cabala.  Scan  through  this  list  and  compare  your  program  to  each  of 
the  questions  asked  here. 

6.1  Data  Reference  Errors 

1 .  Is  a  variable  referenced  whose  value  is  unset  or  uninitialized? 

2.  Are  there  any  variables  which  are  not  referenced?  Is  this  because  you 
misspelled  a  variable  name? 

3*  Are  all  array  references  within  the  bounds  of  the  array  size? 

4.  Are  all  array  references  selected  with  an  integer  subscript?  Some 
languages  allow  REAL  type  variables  as  array  subscripts  (it  is  best  to 
avoid  this  'feature').  If  you  do  this,  is  the  subscript  what  you  expect  it 
to  be? 

5.  Is  a  variable  being  assigned  a  value  which  does  not  match  its  type? 

6.  Are  there  any  'off-by-one'  errors  in  referencing  array  elements? 

7.  If  an  array  is  referenced  in  several  procedures  or  subroutines,  is  it 
defined  identically  in  all  sections? 

6.2  Data  Declaration  Errors 

1.  Have  all  variables  been  explicitly  declared?  Undeclared  arrays  can  be 

misinterpreted  by  some  compilers  as  functions. 

2.  Is  each  variable  declared  to  have  the  correct  type? 

3*  Are  all  variables  and  arrays  initialized  properly?  I.e.,  do  the  values 
assigned  to  each  variable  agree  with  the  types  of  each  variable? 

4.  Are  there  any  variables  with  similar  names  (e.g.,  POINT  and  POINTS)?  This 
is  not  necessarily  an  error,  but  it  is  a  sign  that  names  may  have  been  con¬ 
fused  or  misspelled  somewhere  in  the  program. 

6.3  Computation  Errors 

1.  Are  computations  done  with  variables  having  inconsistent  types  (e.g.,  mul¬ 
tiplication  of  variables  of  type  'character')? 

2.  Are  there  any  mixed-mode  computations  (INTEGER  and  REAL)?  This  Is  not 
necessarily  an  error,  but  is  the  computation  result  of  the  expected  type9 


WsrTf 


THE  PROGRAM  CHECKLIST 


Computation  Errors 


3. 

4. 

5. 

6. 

7. 

8. 

9. 

6.4 
1. 

2. 

3. 

4. 

5. 

6.5 
1. 


Do  computations  with  arrays  have  the  proper  matching  lengths  where 
required?  E.g.,  a  matrix  product  requires  a  match  of  the  'inner  dimension' 
of  the  arrays. 

Are  any  divisions  performed  with  divisors  very  close  to  zero?  Does  this 
affect  the  validity  of  the  computations? 

Is  it  possible  for  a  divisor  to  be  zero?  Is  it  possible  that  some  func¬ 
tions  employ  a  zero  divisor  (e.g.,  the  MOD  function,  the  arctangent  func¬ 
tion,  etc.)? 

Are  there  any  consequences  of  the  fact  that  digital  computers  rarely 
represent  decimal  numbers  exactly?  I.e.,  1/3  +  1/3  +  1/3  does  not  equal' 
1.0. 


Can  a  variable  go  outside  its  meaningful  range?  For  example,  can  a  vari¬ 
able  measuring  probability  ever  be  greater  than  1.0? 

Is  the  assumption  of  the  order  of  evaluation  correct  for  expressions  which 
contain  more  than  one  operator? 


Are  there  any  invalid  uses  of  integer  arithmetic?  For  instance,  if  I  is  an 
integer  variable,  2*1/2  is  equal  to  I  only  if  I  is  even  and  only  if  the 
multiplication  is  done  first. 


Comparison  Errors 


Are  there  any  comparisons  of  variables  of  incompatible  types?  I.e.,  are 
characters  compared  with  reals? 

Are  there  any  mixed  mode  comparisons? 

Are  the  comparison  operators  correct?  Most  of  the  difficulty  arises  in  the 
combined  use  of  'and',  'or',  and  'not.' 

Are  the  operators  of  logical  expressions  of  type  logical?  For  example,  to 
determine  if  I  is  between  the  values  2  and  10,  the  correct  expression  is 
(2<I)A(I<10),  not  (2<I<10). 

Are  there  any  comparisons  of  numbers  with  fractional  parts  in  which  trunca¬ 
tion  errors  play  a  role? 

Construct  Errors 

Is  it  possible  that  certain  entry  conditions  will  prevent  execution  of  your 
program?  For  example,  in  the  loop 

while  "not  found" 

S. 

what  happens  if  "found"  is  initially  true? 
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2.  Can  the  number  of  selections  in  a  Multiple  Choice  construct  ever  exceed  the 
number  of  possibilities  you've  allowed  for? 

3.  Will  every  loop  eventually  terminate? 

4.  What  are  the  consequences  of  an  Indeterminate  Iteration  going  all  the  way 
to  the  bitter  end?  For  example,  in  the  following  program  fragment  (a  loop 
controlled  by  a  compound  Boolean  expression),  what  happens  if  "found"  never 
becomes  true? 

i  becomes  1 . 

while  (  (i  <  tablesize)  and  (not  found )){ 

S. 

i  becomes  i  +  1 . 

} 

5.  Are  there  any  "off-by-one"  errors  (too  many  or  too  few  iterations)? 

6.  Is  there  a  corresponding  ending  bracket  to  every  opening  bracket  (if  the 
programming  language  uses  them)?  Are  the  program  statements  grouped  prop¬ 
erly  with  brackets? 

7.  Are  there  any  non-exhaustive  decisions?  E.g.,  if  a  variable  is  supposed  to 
have  one  of  the  values  1,  2,  or  3,  do  you  assume  that  the  value  must  be  3 
if  it  is  not  1  or  2?  Is  this  valid,  particularly  for  program  input? 

6.6  Interface  Errors 

Some  programming  languages  allow  programs  to  be  broken  into  smaller,  more 
manageable  parts  called  subroutines  or  modules.  The  interface  between  modules 
can  sometimes  be  a  source  of  errors. 

1 .  Is  the  number  of  arguments  received  by  a  module  the  same  as  the  number 
sent?  Are  they  in  the  correct  order? 

2.  Are  there  any  unused  arguments? 

3.  Do  the  attributes  of  each  passed  variable  match  the  attributes  of  the 
received  variable?  For  example,  is  a  simple  REAL  variable  passed  to  a 
module  which  expects  an  array? 

4.  Are  the  units  of  each  passed  variable  correct?  For  example,  is  the  passed 
value  expressed  in  degrees  while  the  expected  value  is  expressed  in  radi¬ 
ans? 

5.  Does  the  number  of  variables  passed  by  a  module  equal  the  nunber  of  vari¬ 
ables  expected  by  the  called  module?  Are  they  in  the  correct  order? 

6.  Are  the  number,  order,  and  type  of  variable  correctly  passed  to  'built-in' 
functions?  For  example,  does  the  arctangent  function  require  one  or  two 
arguments?  If  two,  which  one  is  the  divisor? 


& 
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Interface  Errors 


7.  Does  a  subroutine  alter  a  value  which  is  supposed  to  be  only  input?  For 
example,  is  an  argument  defining  the  length  of  an  array  used  as  the  index 
to  a  determinate  loop?  What  is  it's  value  on  return  from  the  module? 

8.  Are  global  variables  referenced  the  same  way  in  all  modules?  In  FORTRAN, 
for  example,  are  the  variables  listed  in  COMMON  blocks  the  same  everywhere? 

9.  Are  constants  ever  passed  to  subprograms?  For  example,  the  FORTRAN  state¬ 
ment 

CALL  SUB(A, 3) 

can  be  dangerous,  because  if  SUB  assigns  a  value  to  the  second  argument,  3 
will  no  longer  be  3!  (This  has  given  rise  to  the  old  saw:  'All  constants 
are  variable.') 

10.  Will  the  program,  subroutine  or  module  eventually  terminate? 

6.7  Input/Output  Errors 


Many  languages  permit  several  ways  to  access  files.  If  you  are  using  these 
features,  check  your  program  with  this  list: 

1 .  Are  file  attributes  correctly  declared? 


2.  Are  the  attributes  on  the  OPEN  statement  correct? 

3.  Does  the  format  specification  agree  with  the  READ  statement? 

4.  Are  arrays  declared  to  be  large  enough  to  contain  all  the  information  to  be 
read? 

5.  Have  all  files  been  opened  before  use? 

6.  Are  end-of-file  conditions  detected  and  handled  correctly? 

7.  Are  Input/Output  error  conditions  handled  correctly? 

8.  Are  there  spelling  or  grammar  errors  in  any  messages  written  by  the  pro¬ 
gram?  Are  the  messages  intelligible? 

6.8  Miscellany 

1.  If  the  compiler  you  are  using  has  a  'post  mortem  dump'  switch,  set  it  to 
'on.'  This  switch  will  allow  the  computer  to  print,  at  the  time  of  your 
program's  death,  the  values  of  your  variables  according  to  their  types 
(rather  than  in  their  octal  or  hexadecimal  representation).  I.e.,  the 
values  of  REAL  variables  are  printed  in  engineering  notation,  INTEGERS  are 
printed  as  Integers,  LOGICALs  are  printed  as  TRUE  or  FALSE,  etc.  The  clues 
offered  by  this  tool  are  extremely  helpful.  (Note  that  the  procedures  of 
Chapter  4  of  [Roth83]  have  this  switch  set  'on'.) 
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2.  Does  the  compiler  cross  reference  map  Indicate  variables  which  are  unused 
or  referenced  only  once?  This  may  not  be  an  error,  but  might  point  to 
misspellings. 

3.  Are  the  attributes  which  the  compiler  assigns  to  each  variable  the  ones  you 
expected? 

4.  Did  the  compiler  produce  any  'warning'  messages  (assuming  of  course  that 
you  have  a  successful  compile)?  These  messages  point  to  potential  prob¬ 
lems. 

5.  Is  the  program  or  module  sufficiently  'robust?'  That  is,  does  it  cheek  its 
input  for  validity,  or  can  it  be  killed  with  a  'reasonable'  number? 

6.  Does  the  program  NOT  do  something  that  you  expected  it  would? 


7.  DEBUGGING 


Your  program  has  finally  'compiled'  and  *go'ed,  but  produces  erroneous  output. 
This  is  caused  by  'bugs'  which  must  be  found  and  exterminated  in  order  for  the 
program  to  be  'correct.'* 

This  chapter  provides  an  approach  to  bug  eradication;  don't  forget  the  checklist 
in  the  previous  chapter.  Many  of  the  following  suggestions  are  from  [Myer79l. 

7.1  Think 


The  building  size  necessary  to  contain  a  computer  with  the  power  of  your  brain 
would  exceed  that  of  the  Etapire  State  Building.  You  may  not  be  able  to  calcu¬ 
late  quickly,  but  you  can  intuit.  You  should  be  able  to  debug  most  of  your  pro¬ 
grams  without  going  near  a  computer. 

7.2  If  You  Reach  An  Impasse 

-  sleep  on  it.  Your  subconscious  mind  has  great  potential  for  working  things 
out  while  you're  doing  something  else.  This  is  not  an  excuse  to  catch  a 
few  z's  on  the  job. 

-  describe  the  problem  to  someone  else.  By  making  an  effort  to  tell  a  good 
listener  what  your  program  is  doing,  you  may  discover  the  problem  yourself. 
The  listener  need  not  be  an  expert  on  what  you're  doing,  either;  an 
ignorant  person  can  often  see  the  naked  emperor. 

7.3  Where  There  Is  One  Bug 

there  are  likely  to  be  more.  Examine  the  immediate  vicinity  of  a  bug  for  other 
errors,  since  bugs  usually  result  from  a  misunderstanding  of  what  the  program  is 
supposed  to  do.  Note  also  that  small  pieces  of  code  may  be  veritable  'roach 
hotels;'  there  are  such  things  as  error-prone  modules. 

7.*»  Fix  The  Error 


Don't  just  fix  the  symptom  of  the  error.  Make  sure  that  all  occurrences  of  the 
error  are  fixed,  not  just  this  'just-discovered'  one. 

7.5  Bebugging 

Recognize  that  fixing  a  program  is  likely  to  introduce  new  bugs,  because  the 
entire  concept  of  the  program  may  not  be  fresh  in  your  mind. 

7.6  Faulty  Bug  Repair 

Your  'fix'  may  be  wrong!  The  probability  of  a  correct  fix  decreases  with  the 
size  of  the  program,  and  varies  with  the  size  of  the  fix.  Bug  repair  should 
place  you  mentally  in  the  program  design  stage. 

*  If  you  would  rather  wave  your  hand  and  say,  'Yeah,  I  know  it  doesn't  work 
for  such-and-such  a  case,  but  so  what!?'  then  the  bug  is  known  as  a 
' feature. ' 


Use  Debugging  Tools 

7.7  Use  Debugging  Tools 


DEBUGGING 


only  as  a  second  resort.  An  'interactive  debugger'  requires  that  you  learn 
another  language,  and  only  gives  a  static  picture  of  what  a  program  is  doing. 
Use  them  as  an  adjunct  to,  but  not  a  replacement  for,  thought. 

7.8  Avoid  Experimentation 

Experiment  only  as  a  last  resort.  A  common  mistake  made  by  novices  is  to 
attempt  to  solve  the  problem  by  changing  the  program.  For  example,  nHmmmm...T 
don't  know  what's  wrong,  so  let's  change  this  DO  loop  and  see  what  happens." 
This  kind  of  behavior  has  little  chance  of  addressing  the  actual  problem,  and 
muddies  the  waters  by  introducing  other  errors. 

7.9  Use  an  Octal  or  Hexadecimal  Dump 
never. 
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8.  THE  LANGUAGES 


The  second  most  important  requirement  of  computer  programs  is  that  they  be  legi¬ 
ble  to  humans.*  We  use  programming  languages  when  we  talk  with  computers  because 
English  is  imprecise  (although  we  humans  can  understand  it  even  when  it  is 
misused!).  While  we  maintain  precision  in  our  programming,  our  computer  talk 
must  also  be  understandable  by  hunans,  because  we  are  the  race  which  must  'main¬ 
tain'  the  programs.  This  understandability  is  controlled  first,  by  the  language 
itself,  and  second,  by  the  way  in  which  we  use  the  language. 

The  languages  are  presented  in  the  order  of  decreasing  readability.  Thus,  if 
you  wish  your  programs  to  be  readable,  you  should  rather  write  in  ratfor,  Pas¬ 
cal,  or  C,  and  avoid  BASIC  and  programmable  calculators.  My  viewpoint  is:  the 
language  defines  the  machine.  Hence,  I  may  write  programs  for  a  'Pascal  com¬ 
puter  or  a  'FORTRAN  Engine'  and  not  concern  myself  with  what  the  machine  really 
is. 

An  aside  re  'standard'  languages:  there  are  no  such  things!  Although  a  document 
may  exist  which  describes  the  'standard'  for  a  particular  language,  each  imple¬ 
mentor  includes  the  bells  and  whistles  which  make  his  own  hardware  hum.  So 
there  are  actually  several  'versions'  of  a  standard  language;  moving  from  one 
machine  to  another  can  still  be  a  difficult  proposition. 

8. 1  Languages  Included 

ratfor  Invented  by  Kernighan  and  Plauger  [Kern761  to  provide  reasonable  con¬ 
structs  for  FORTRAN,  rational  fortran  is  really  a  preprocessor  to  the 
FORTRAN  language.  It  allows  all  the  features  of  FORTRAN  to  be 
accessed,  so  the  advantages  of  portability,  universality,  and  relative 
efficiency  of  FORTRAN  are  retained.  In  addition,  it  has  a  'macro' 
expansion  capability,  is  free-format,  and  allows  comments  on  the  same 
line  as  program  text.  This  language  is  easily  read  by  programmers 
with  a  FORTRAN  background. 

ratfor  is  maintained  in  two  versions:  a  source  version  in  ratfor 
itself  (to  permit  maintenance  in  a  reasonable  language)  and  a  source 
version  in  FORTRAN66  to  allow  easy  portability  amongst  computers. 

Pascal  A  simple  language  invented  by  Niklaus  Wirth  [Jens741  specifically  to 
teach  programming,  it  has  found  wide  use  on  microcomputers  as  well  as 
on  mainframes.  Its  structure  'enforces'  good  programming  practices, 
and  catches  many  errors  early  in  the  programming  process,  saving 
debugging  time  and  cost.  It  is  free-format,  and  allows  extreme  flexi¬ 
bility  in  data  structuring,  commenting,  and  recursion  in  both  data  and 
procedures.  The  best  texts  for  learning  Pascal  are  [Wirt73l, 
[Jens74],  [Atki80],  and  [Coop82]. 

This  is  the  language  in  which  the  UNIX  operating  system  is  written. 
It  is  a  free-format,  terse,  language  whose  form  contributed  much  to 
the  style  of  ratfor.  Several  significant  scientific  programs  have 
been  written  in  C.  The  standard  C  manual  is  the  reference  [Kern781. 


*  Readability  supersedes  all  program  requirements  except  correctness 
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8.2 

K- 

•V* 

Ada* 

Bfl 

m 

Algol 

*  * 

APL 

FORTRAN77  The  most  recent  revision  of  FORTRAN.  The  major  updates  it  provides  to 
FORTRAN66  are:  the  CHARACTER  type,  string  processing,  the  block  IF- 
THEN-ELSE,  and  greater  explicit  control  over  Input/Output  (I/O).  Syn¬ 
tax  is  still  somewhat  inconsistent,*  and  the  meaning  of  one  of  the  FOR- 
TRAN66  constructs  has  been  changed(l). 

FORTRAN66  The  'standard'  scientific  language  since  1954.  It  is  the  lingua 
franca  of  the  engineering  community.  All  of  the  'large'  finite  ele¬ 
ment  programs  have  been  written  in  FORTRAN66,  and  almost  all  computer 
vendors  offer  a  FORTRAN66  compiler.  Syntax  of  this  language  is 
unusual,  and  it  forces  you  to  write  your  program  backwards  from  the 
way  you  would  express  it  in  English.  But  'everyone'  knows  it,  so  its 
anomalies  are  generally  overlooked  or  unrecognized.  It  is  to  your, 
advantage  to  write  in  some  other  language  if  possible;  'ratfor'  is  a 
r easonabl e  alternative . 

BASIC  The  most  widespread  microcomputer  language  available  today.  Invented 
to  teach  computing  to  Dartmouth  students,  this  language  is  feature- 
poor,  and  also  forces  you  to  program  backwards  from  your  thinking 
processes. 

8.2  Languages  Excluded 

Ada*  This  language  is  now  under  development.  It  is  the  new  'standard'  DoD 

language  for  embedded  systems  (systems  which  are  part  of  some  other 
machine,  such  as  a  tank,  a  ship,  a  torpedo,  etc.).  Ada  is  an 
'enhanced  derivative'  of  Pascal. 

Algol  This  family  of  languages  preceded  the  development  of  Pascal.  It  is 

the  native  language  of  Burroughs  machines. 

APL  A  very  powerful,  very  cryptic,  language.  The  information  content  of 

APL  programs  is  extremely  dense.  Per  [KellSl],  APL  is  "so  compacted 
that  the  source  code  can  be  freely  disseminated  without  revealing  the 
programmer's  intent  or  Jeopardizing  proprietary  rights."  It  is  impos¬ 
sible  to  maintain  programs  written  in  APL;  one  throws  them  away  and 
redoes  them. 


Assembly 


Assembly  language  is  a  mnemonic  method  for  addressing  the  computer's 
'instruction  set.'  Usually,  it  is  a  very  primitive  language  (note, 
however,  that  the  Burroughs  machine  instructions  are  Algol,  a  'high 
level'  language).  Programming  in  assembly  language  is  undesirable 
because  it  takes  too  much  time  to  do  it,  and  because  it  is  nigh  impos¬ 
sible  to  move  it  to  another  vendor's  machine.  Certainly,  some  pieces 
of  some  programs  should  be  written  in  assembly  language  (data  acquisi¬ 
tion  programs  and  arcade  games,  for  example),  but  because  it  can  take 
ten  times  more  human  time  to  write  assembly  than  some  other  higher- 
order  language,  it  is  rare  that  one  can  recover  an  investment  in 
assembly  language  code  development.  The  code  produced  by  optimizing 
compilers  is  as  good  as,  and  in  many  cases,  better  than,  the  code  one 
could  produce  'by  hand'  using  assembly  language.  Therefore,  assembly 


Ada  is  a  trademark  of  the  Department  of  Defense. 
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language  programming  is  to  be  eschewed. 

COBOL  A  language  designed  primarily  for  business  purposes,  and  not  very  use¬ 
ful  for  engineering  purposes. 

FORTH  A  'threaded'  language  in  some  use  on  microcomputers,  which  forces  the 
programmer  to  act  as  part  of  the  compiler.  This  language  immortalizes 
reverse  Polish  notation  (alias  'Okie  code')  in  programs;  FORTH  pro¬ 
grams  are  hard  to  read. 

Lisp  The  de  facto  language  of  the  artificial  intelligence  community.  It  is 

not  designed  for,  and  hence  has  seen  little  use  in,  number  crunching. 

Modula-2  Niklaus  Wirth's  successor  to  Pascal.  It  removes  many  of  the  difficul¬ 
ties  associated  with  program  development  in  Pascal,  but  distribution 
of  this  language  has  not  been  wide-spread  to  date. 

PL/I  IBM's  answer  to  the  FORTRAN/COBOL  dichotomy.  An  enormous  and  complex 

language.  Never  really  caught  on  in  the  engineering  community. 

TI59  The  Texas  Instruments'  model  59  Programmable  calculator  (as  well  as 

other  models  by  Texas  Instruments,  Hewlett-Packard,  Sharp,  etc.) 
allows  programs  to  be  written  in  terms  of  keystrokes.  The  languages 
are  very  close  to  assembly  language;  hence,  the  remarks  directed  at 
assembly  languages  are  also  applicable  here.  Although  I  planned  to 
include  examples  of  the  six  constructs  expressed  in  this  language, 
they  were  so  unreadable  that  I  omitted  them.  Because  of  the  recent 
advances  in  miniaturization  of  electronic  components,  shirt-pocket 
computers  are  now  available  (with  BASIC  interpreters)  for  the  price  of 
a  calculator.  Programmable  calculators  are  antiques  (which  is  not  to 
say  that  they  are  not  useful!). 

others  Languages  which  are  either  dead,  distasteful,  experimental,  unknown  to 
me,  or  little  used  for  engineering  work.  This  includes  languages  like 
MAD,  CLU,  Gypsy,  SNOBOL,  Euclid,  Simula,  RPG,  Simscript,  etc.,  etc. 
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9.  RATFOR 


Blocks  of  statements  In  ratfor  are  grouped  with  braces  '  { '  and  '}'.  On  CDC  sys¬ 
tems,  these  may  be  entered  as  either  braces  or  brackets  '[  ]'.  Unfortunately, 
the  CDC  character  set  is  not  common  across  all  the  printers,  so  brackets  may 
appear  as  either  brackets  or  braces. 

ratfor  is  a  free-format  language.  It  is  a  good  idea  to  indent  your  code  so  that 
the  logical  structure  of  the  program  is  reflected  by  the  layout  of  the  text  on  a 
page. 

There  are  three  ways  to  provide  commentary  in  ratfor  programs:  with  a  *C'  or 
in  column  1,  or  with  a  anywhere.  You  should  start  your  program  text  in 
column  2  or  3  so  that  a  'C'  or  a  '*'  in  column  1  is  not  interpreted  as  a  comment 
line.  Blank  lines  are  ignored  by  ratfor,  so  you  may  also  use  these  to  provide 
additional  'white  space.' 

9.1  Variable  Names 

Names  of  variables  are  limited  by  the  limitations  of  the  target  FORTRAN  (since 
ratfor  preprocesses  text  to  FORTRAN).  Normally,  this  is  a  maximum  of  six  char¬ 
acters.  Every  variable  should  be  declared  to  be  of  a  specific  type,  and  its 
purpose  defined.  This  may  be  done  easily  through  the  type  definitions  REAL, 
INTEGER,  LOGICAL,  CHARACTER,  COMPLEX,  and  DOUBLE,  and  the  on-the-same-line  com¬ 
mentary  ratfor  provides.  It  is  inadvisable  to  rely  on  FORTRAN  to  automatically 
type  your  variables  for  you. 

The  type  definitions  may  also  be  used  to  declare  the  dimensions  of  arrays,  so 
the  DIMENSION  statement  is  not  needed  in  ratfor  or  FORTRAN  programs. 

ratfor  uses  the  'reserved  word'  concept?  i.e.,  there  are  several  words  which  are 
special  to  ratfor.  Thus,  in  ratfor  programs,  you  may  not  use  any  of  the  follow¬ 
ing  words  for  variable  names:  break,  do,  else,  if,  next,  repeat,  until,  and 
while. 
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9*2  Boolean  Expression 


X  and  Y  are  assumed  to  be  variables  of  the  same  type. 


X  is  greater  than  Y: 


X  >  Y 


X  is  greater  than  or  equal  to  Y: 


X  >=  Y  or  X  =>  Y 


X  is  equal  to  Y: 


X  ==  Y 


X  is  not  equal  to  Y: 

X  *=  Y  (on  CDC) 

X  !=  Y  (on  UNIX) 


X  is  less  than  Y: 


X  <  Y 


X  is  less  than  or  equal  to  Ys 


X  <s  Y  or  X  s<  Y 


9 . 3  Assignment  Statement 


X  =  expression 
(Same  as  FORTRAN.) 
9.4  Conditional 


if  (B)  { 

51 

52 

}  #  end  if 


Alternative 


RATFOR 


9.5  Alternative 

if  (B)  { 

51 

52 

} 

else  { 

53 

54 

}  #  end  if 

9.6  Multiple  Choice 

The  clearest  way  to  select  from  multiple  possibilities  is  with  a  series  of  If- 
Then-Else's. 

if  (  B1  ) 

51 

else  if  (  B2  ) 

52 

else  if  (  B3  ) 

53 

•  •  • 
else 
Sn 

9.7  While 

while  (B)  { 

51 

52 

}  #  end  while 

9.8  Repeat-Untll 

repeat  {  #  until  (B) 

51 

52 

}  until  (B) 


where  it  is  helpful  to  a  reader  to  have  the  termination  condition  of  the  loop 
presented  in  commentary  on  the  Repeat  line. 
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9.9  Determinate  Iteration 


The  ratfor  'Do'  is  identical  to  the  FORTRAN  'Do',  with  the  exception  that  ratfor 
requires  no  label*  as  a  loop  terminator.  If  the  target  language  is  FORTRAN66, 
then  the  Do  index  may  only  increase;  if  the  target  is  FORTRAN77,  the  index  may 
increase  or  decrease.  (Hence,  a  decreasing  index  is  more  clearly  handled  with  a 
While  or  a  Repeat-Until.) 

do  I  =  J,K,M  { 

51 

52 

}  #  end  do 
9.10  Altered  Loops 

9.10.1  Loop  Exit 

Although  any  loop  (While,  Repeat-Until,  Determinate)  may  be  'exit'ed,  we  illus¬ 
trate  with  a  While. 

while  (B1 )  { 

51 

if  (B2) 
break 

52 

}  #  end  while 

9.10.2  Loop  Redo 

Although  any  loop  (While,  Repeat-Until,  Determinate)  may  be  ' re-do 'ed,  we  illus¬ 
trate  with  a  While. 

while  (B1 )  { 

51 

if  (B2) 
next 

52 

}  #  end  while 


*  In  general,  the  only  labels  you  need  in  a  ratfor  program  are  those  on  FORMAT 
statements. 
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Generalized  Iteration 
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9.11  Generalized  Iteration 

ratfor  provides  a  generalized  iteration  scheme  called  the  For.  It  has  the  form 

for  (  initialize;  B;  reinitialize  )  { 

51 

52 


which  is  the  equivalent  of 

initialize 
while  (B)  { 

51 

52 

reinitialize 

1 

The  For  may  be  easier  to  understand  in  some  circumstances  because  it  keeps  all 
terms  which  control  the  loop  on  a  single  'line,' 
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9.12  rat for  Summary 

Conditional: 

Alternative: 

Multiple  Choice: 

While: 

Repeat-Until: 

Determinate: 

Generalized  loop: 

Loop  Exit: 


if  (B) 

S 

if  (B1 ) 

51 
else 

52 

if  31 

51 

else  if  B2 

52 

else  if  B3 

53 

else  if  84 

54 

while  (B) 

• *  » 

S 

repeat 

S 

until  (B) 

do  i  s  j,k,m 
S 

for  (  init;  Bj  re-init  ) 
S 

break 


Loop  Redo 


next 


Program  Header 


RATFOR 


The  ratfor  program  header  requires  a  comment  line  whioh  identifies  the  program, 
followed  by  the  seme  information  as  that  required  in  a  FORTRAN  program.  On  CDC 
machines,  this  has  the  form 


#  name  -  a  line  to  demonstrate  the  header  format 
program  name(  filel,  file2,  ...  ) 


declaration  of  variable  types 


9.14  Running  a  ratfor  Program 

Assume  you  have  composed  your  program  with  a  text  editor  on  a  CDC  computer  and 
the  code  is  in  <p>.  Assume  further  that  you  have  already  executed  the  commands 

ATTACH , ROTH , CCLLIB , ID=CSPR . 

LIBRARY, ROTH. 

(you  must  execute  these  commands  only  once).  Then  the  program  in  <p>  may  be  run 
with  the  command 

rr,p. 

(see  [Roth83l). 
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9.15  Examples 

I  will  assume  that  rat for  is  preprocessing  programs  for  FORTRAN66. 

9.15.1  Temperature  Conversion 

#  ftoc  -  fahrenheit  to  Celsius  conversion,  -40f  to  lOOf 
program  main (output) 

real  f,c  #  fahrenheit,  Celsius  temperatures 
f  =  -40.0 

whileC  f  <=  100.0  )  { 
c  =  5.0*(f  -  32.0)/9.0 
print  *,  f,c 


Examples 


RATFOR 


9*15*2  Nearest  Points  Assume  that  the  three  coordinates  of  each  point  are 
typed  on  a  single  line  of  <input>. 

#  near  -  find  2  of  10  points  which  are  closest  neighbors 
program  near(  input,  output) 
real  x(10,3)  #  array  of  points 

real  oldd,newd  #  distances  between  points 
integer  l,j,k,m,n  #  counters 

read  •,((x(i,j),J=1,3),i=1,10)  #  free  format  read 

#  first  distance t  from  1  to  2 

oldd  s  sqrt(  (x(1,1)-x(2,1))*«2 
+  (x(1 ,2)-x(2,2))**2 
+  (x( 1 ,3)-x(2,3) )**2  ) 

do  i  =  1,  9  { 

ipl  s  i  +  1 

do  j  =  ipl,  10  { 

newd  =  sqrt(  (x(i,1)-x(J,1))»»2 
+  (x(i,2)-x( j,2))**2 
+  (x(i,3)-x( j,3))**2  ) 

if(  newd  <=  oldd  )  {  #  found  a  new  nearest  neighbor 

m  s  i 

n  s  J 

oldd  s  newd 
} 

)  #  end  do  J 

}  #  end  do  i 

print  closest  points  are  at  ' 
print  •,  (x(m,i) ,1=1 ,3) ,  (x(n,i),i=1, 3) 
end 
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Examples 


9.15.3  Count  the  'A's  Approach:  assume  that  each  line  has  at  most  132  Qharao- 
ters,  and  read  an  entire  line  at  a  time.  Also  assume  that  character  variables 
can  be  represented  as  integers  (this  works  on  CDC  machines).  If  the  line  is 
shorter  than  132  characters,  FORTRAN  will  act  as  if  the  nonexistent  characters 
were  blanks.  If  the  line  is  longer  than  132  characters,  the  excess  will  be 
(silently)  truncated  (perhaps  not  the  best  solution), 

#  counta  -  count  the  number  of  times  'A'  appears  in  <input> 
program  counta(  input, output, tape5=input  ) 
integer  count,  line (132),  i 

integer  eof  #  end-of-file  funotion  (cdc  supplied) 
count  r  0 

read (5,1)  (line(i) ,i=1 , 132) 

1  format (1 32a 1) 
while (  eof (5)  ==  0  )  { 
do  i  =  1,  132 

if (  line(i)  ==  'A'  ) 
count  =  count  +  1 
read (5, 1 )  (line(i) ,i=1 , 132) 

} 

print  •, count 


10.  PASCAL 


Blocks  of*  statements  in  Pascal  are  grouped  with  the  symbols  'begin'  and  'end'; 
you  might  think  of  them  as  'fat  brackets.' 

Pascal  is  a  free-format  language.  That  is,  you  may  begin  your  code  in  any 
column.  It  is  a  good  idea  to  indent  your  code  so  that  the  logical  structure  of 
the  program  is  reflected  by  the  layout  of  the  text  on  a  page. 

The  symbol  '(*'  or  '{'  opens  a  comment,  and  the  next  occurrence  of  the  symbol 
'*)'  or  closes  a  comment.*  Comments  may  occur  anywhere  and  may  extend  across 
line  boundaries.  Blank  lines  are  ignored  by  Pascal,  so  you  may  also  use  these 
to  provide  additional  'white  space.' 

10.1  Variable  Names 

As  in  English,  Pascal  does  not  limit  the  number  of  characters  in  a  variable 
name.  Most  implementations  of  the  language  only  guarantee  that  the  first  8 
characters  are  significant.  That  is, 


a  fairly 1 ongname 

-'"'afairlylongnamewithstuf  fontheer.J 

probably  refer  to  the  same  variable.  Although  it  is  possible  to  use  any  number 
of  characters,  limit  yourself  to  at  most  12. 


Note  that  Pascal  forces  you  to  declare  the  types  of  your  variables  before  you 
use  them.  So  while  you  are  declaring  their  types,  you  might  just  as  well  Insert 
some  commentary  which  declares  their  purposes. 


Pascal  uses  the  'reserved  word'  concept;  i.e.,  many  words  are  special  to  Pascal. 
Thus,  in  Pascal  programs,  you  may  not  name  a  variable  'if,  while,  repeat, 

etc.  A  full  list  of  the  Pascal  reserved  words  is  given  in  tJens741. 


*  The  braces  are  not  available  on  CDC. 
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10.2  Boolean  Expression 

X  and  Y  are  variables  of  the  sane  type. 

X  is  greater  than  Y:  -- 

X  >  Y 

X  is  greater  than  or  equal  to  Ys 
X  >=  Y 

X  is  equal  to  Ys 
X  =  Y 

X  is  not  equal  to  Ys 
X  <>  Y 

X  is  less  than  Ys 
X  <  Y 

X  is  less  than  or  equal  to  Ys 
X  <r  Y 

10.3  Assignment  Statement 
X  s=  expression 

10.4  Conditional 
if  B  then  begin 

si; 

S2 

end; 

The  semicolon  is  used  in  Pascal  to  .join  statements.  'begin'  and  'end'  are 

reserved  words'  which  serve  the  same  function  as  brackets  in  ratfor. 


V  ».  _ 


i *. 


Alternative 
40r5  Alternative 
if  B  then  begin 


PASCAL 


A. '  »  " 

SI; 

v  *  • 

&:■ 

S2 

end 

5H 

else  begin 

UWm 

S3; 

S4 

S-:-v; 

end  (•if*) 

m 

10.6  Multiple  C 

There  are  two  ways  of  selecting  from  multiple  alternatives:  the  case  and  the 
if  ...  el3e  if  constructs. 

10.6.1  The  Case 

case  I  of 
LI:  SI; 

L2:  S2; 

L3:  S3; 

•  •  • 

Ln:  Sn 

end  (*case»); 

The  labels  'Ll',  'L2',  'L3',  •••,  'Ln'  are  the  legitimate  values  which  'I'  can 
assume. 

10.6.2  if  ...  else  if 

if  B1  then 

51 

else  if  B2  then 

52 

else  if  B3  then 

53 

else  if  B4  then 
S 4 ( *  end  if  »); 
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10.7  While 

while  B  do  begin 
SI? 

S2 

end (•while*); 


Repeat-Until 


repeat  (•  until  B  •) 
SI? 

S2 

until  B; 


This  construct  is  made  clearer  to  a  reader  by  putting  the  termination  condition 
in  commentary  on  the  Repeat  line, 


10.9  Determinate  Iteration 
10.9.1  Increasing  Index 


for  I  :=  J  to  K  do  begin 

si; 

S2 

end(»for») ; 


10.9.2 


Index 


for  I  :=  J  downto  K  do  begin 

si? 

S2 

end (•for1) ; 


In  these  constructs,  J  is  the  initial  value  of  I,  and  K  is  the  terminal  value. 
Note  that,  if  the  initial  condition  satisfies  the  test  for  termination,  the  loop 
will  not  be  executed.  In  Pascal,  the  increment  to  a  loop  may  be  +1  or  -1.  No 
other  values  are  possible.  For  increments  other  than  +1  or  -1,  use  a  While  or 
Repeat-Until  construct. 
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10.10  Alterec.  Loops 


10.10.1  Loop  Exit  Although  any  loop  (While,  Repeat-Until ,  Determinate)  may  be 
'exit'ed,  we  illustrate  with  a  While. 


label  13; 


while  B1  do  begin 

Cl  • 

>-  1 » 

if  B2  then  goto  13? 

S2 

end(*while*) ; 

13:  ( "continue  program* ) 

Although  this  construct  is  possible,  and  indeed  works,  most  Paseal-ers  would  not 
recommend  it.  Clearer,  and  more  in  keeping  with  the  style  of  the  language  would 
be  something  like 

keepon  :=  true;  (•  comment  explaining  'keepon'  *) 
while  B1  and  keepon  do  begin 
SI; 

if  B2  then 

keepon  :=  false 
else  begin 
S2 
end 

end(*while*) ; 

because  it  presents  all  the  loop  terminators  in  a  single  statement  (the  While 
statement)  and  it  presents  the  logic  in  explicitly  logical  terms  rather  than  in 
terms  of  transportation  ('goto  label'). 
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Altered  Loops 


10.10.2  Loop  Redo  Although  any  loop  (While,  Repeat-Hntil, 
*  re-do 'ed,  we  illustrate  with  a  While, 

label  13; 


while  B1  do  begin 
SI; 

if  B2  then  goto  13; 

S2 

13:  end(*while*); 

Although  this  construct  is  also  possible,  and  indeed  works, 
it.  Rather,  the  following  is  preferred* 

while  B1  do  begin 
31; 

if  B2  then  begin 

(*  empty  'begin  end'*  a  do-nothing!  •) 
end 

else  begin 
S2 
end 

end(*while»); 
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10.11  Pascal  Summary 

Conditional: 

Alternative: 

Multiple  Choice: 

While: 

Repeat-Until: 

Determinate: 

Loop  exit: 


if  B  then 
S; 

if  B1  then 
SI 
else 
S2; 

case  I  of 
BH  SI 5 
B2:  S2; 
B3:  S3; 
B4:  S4 
end; 


while  B 
S; 

repeat 
S 

until  B; 

for  I  :s  J  (toldownto)  K  do 
S; 

additional  loop  control  variable 


if  B1  then 

51 

(or)  else  if  B2  then 

52 

else  if  B3  then 

53 

else  if  B4  then 
S4; 
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additional  Conditional  construct 
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10.12  Program  Header 

The  Pascal  program  header  requires  the  following  information: 
program  name(  filel,  file2,  ...  ); 
label  declarations  (if  any) 
constant  declarations  (if  any) 
type  definitions  (if  any) 
variable  definitions  (if  any) 
begin 

The  program  is  closed  with  an  'end'  whioh  matches  the  opening  'begin',  followed 
by  a  period.  That  is,  the  last  symbol  in  the  program  is  'end.'. 

10.13  Running  a  Pascal  Program 

Assume  you  have  composed  your  program  with  a  text  editor  on  a  CDC  computer  and 
the  code  is  in  <p>.  Assume  further  that  you  have  already  executed  the  commands 

ATTACH , ROTH , CCLLIB , IDrCSPR . 

LIBRARY, ROTH. 

(you  must  execute  these  commands  only  once).  Then  the  program  in  <p>  may  be  run 
with  the  commands 


Examples 
10.14  Examples 
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10.14.1  Temperature  Conversion 
program  ftoc ( output ) ; 

(•.  ftoc  -  fahrenheit  to  Celsius  conversion,  -40f  to  lOOf  •) 
var  f,c:  real;  (■  fahrenheit,  Celsius  temperatures  •) 
begin 

f  :=  -40.0; 

while  f  <=  100.0  do  begin 
c  :=  5.0*(f  -  32.0)/9.0; 
writeln(f ,c); 
f  :=  f  +  5.0 
end  (*while*) 


end 


PASCAL 


Examples 


10.14.2  Nearest  Points  Assume  that  the  three  coordinates  of  each  point 
typed  on  a  single  line  of  <input>. 

program  near(  input,  output); 

(•.  near  -  find  2  of  10  points  which  are  closest  neighbors  *) 
var 

xs  arrayfl. .10, 1. .31  of  real;  (•  points  in  spaoe  *) 
oldd,newd:  real;  (•  distances  between  points  •) 

i,j,k,m,n:  integer;  (•  counters  •) 

begin 

for  i  :=  1  to  10  do 
for  j  :=  1  to  3  do 
read(x[i, j]); 

( *  first  distanoet  fro*  1  to  2  *) 
oldd  :=  sqrt(  sqr(x[ 1 , 1]-x[2, 1]) 

+  sqr(x[1,2]-x[2,2]) 

+  sqr(x[1,3l-x[2,3l)  ); 

for  i  :=  1  to  9  do 

for  j  :=  i+1  to  10  do  begin 

newd  ;=  sqrt(  sqr(x[i, 1 ]-x[J, 11) 

+  sqr(x[i,2]-x[ j,2]) 

+  sqr(xCi,3l-x[j,3l)  ); 

if  newd  <=  oldd  then  begin  (*  found  a  new  nearest  neighbor*) 
m  :=  i; 
n  s=  j; 
oldd  :=  newd 
end  (•  if  *) 
end  (*  for  j  *); 

write  ('  closest  points  are  at  '); 
for  i  ;=  1  to  3  write(x[m,i]) ; 
for  i  :=  1  to  3  write(x[n,il); 
writeln 
end. 


are 
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10.14.3  Count  the  'A's 

P 

program  counta(  input, output  ) 

»  * 

& 

var  chs  char; 

V.* 

aV 

count:  integer; 

begin 

count  : =  0; 

while  not  eof  do  begin 
read(ch); 
if  ch  =  'A'  then 
count  :=  count  +  1 
end(*  while  •); 
writeln( count) 
end. 

In  Pascal,  'eof'  is  a  Boolean  function  which  tests  for  'end-of-file. ' 


11.  c 

Blocks  of  statements  in  C  are  grouped  with  braces  '  f '  and  '}'.  As  of  this  writ¬ 
ing,  no  C  compiler  exists  for  CDC  computers;  it  is  a  language  which  lives  on 
many  other  brands,  however.  These  include  DEC's  PDP-11/xx  and  VAX  11/7xx,  IBM 
machines  of  various  sizes,  and  many  microcomputers. 

C  is  a  free-format  language.  That  is,  you  may  begin  your  code  in  any  column. 
It  is  a  good  idea  to  indent  your  code  so  that  the  logical  structure  of  the  pro¬ 
gram  is  reflected  by  the  layout  of  the  text  on  a  page.  Blank  lines  are  ignored 
by  C,  so  you  may  also  use  these  to  provide  additional  'white  space.' 

Comments  may  occur  anywhere  and  may  extend  across  line  boundaries.  The  symbol 
'/*'  opens  a  comment,  and  the  next  occurrence  of  the  symbol  '•/'  closes  a  com¬ 
ment. 

11.1  Variable  Names 

C  uses  the  'reserved  word'  concept;  i.e.,  there  are  a  large  number  of  words 
which  are  special  to  C.  Thus,  in  C  programs,  you  may  not  name  a  variable  'if,' 
'while,'  'break,'  etc.  A  full  list  of  the  C  reserved  words  is  given  in 
[Kem78]. 

Only  the  first  8  characters  in  a  name  are  significant.  More  may  be  used,  how¬ 
ever.  A  practical  limit  is  12  characters. 

11.2  Boolean  Expression 

X  and  Y  are  variables  of  the  same  type. 

X  is  greater  than  Yj 
X  >  Y 

X  is  greater  than  or  equal  to  Y: 

X  >=  Y 

X  is  equal  to  Ys 
X  ==  Y 

X  is  not  equal  to  Y: 

X  !=  Y 

X  is  less  than  Ys 
X  <  Y 

X  is  less  than  or  equal  to  Y: 

X  <=  Y 
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Assignment  Statement 


11.3  Assignment  Statement 


X  =  expression  ; 


Note  that  in  C,  statements  are  terminated  with  a  semicolon,  rather  than 
as  in  Pascal  (a  subtle,  but  sometimes  painful,  difference). 


loined 


11.4  Conditional 


if  (B)  { 

Si; 

32; 

}  /•  end  if  */ 


11.5  Alternative 


if  (B)  { 

SI; 

S2; 

) 

else  { 

S3; 

S4* 

}  /•  end  if  •/ 


»  _•  A  .>  A.N  A  ’  .  •  k  *  ■  •  »  ’k  *  ■ 
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Multiple  Choice 


While 


C 
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11.10  Altered  Loops 

11.10.1  Loop  Exit  Although  any  loop  (While,  Repeat-Until,  Determinate)  may  be 
'exit'ed,  we  illustrate  with  a  While. 

while  (B1)  { 

SI; 

If  (B2) 
break; 

S2; 

}  /*  end  while  */ 

11.10.2  Loop  Redo  Although  any  loop  (While,  Repeat-Until ,  Determinate)  may  be 
're-do'ed,  we  illustrate  with  a  While. 

while  (B1 )  { 

SI; 

if  (B2) 
continue 
S2; 
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11.11  C  Summary 


wn 


Conditional : 


Alternative: 


Multiple  Choice: 


While: 


if  (B) 


if  (B) 

SI; 

else 

S2; 

switch  (I)  of  { 
case  B1:  SI; 
case  B2:  S2; 
case  B3:  S3; 
case  B4:  S4; 


while  (B) 
S; 


if  (B1 ) 

SI  • 

else  if  (B2) 
S2; 

else  if  (B3) 
S3* 

else  if  (B4) 
S4; 


Repeat-Until: 


Determinate: 


Loop  exit: 
Loop  redo: 


while  (B); 

for(init;  B;  reinit) 
S; 

break 

continue 


pro 


I 


c 
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11.12  Program  Header 


C  is  not  yet  available  on  CDC  machines, 
is 


On  a  UNIX  system,  the  C  program  header 


main  () 

{ 

declaration  of  variable  types 


Note  that  the  program  must  be  closed  with  a  right  brace  ' } '  to  match  the  opening 
left  brace  of  the  program  header. 


11.13  Running  a  C  Program 

Assume  you  have  composed  your  program  with  a  text  editor  on  a  UNIX  system  and 
the  code  is  in  <p.c>.  Then  the  program  in  <p.c>  may  be  run  with  the  commands 


cc  p.c 
a.  out 
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11.14  Examples 

I  assume  in  the  following  examples  that  the  function  ' printf '  is  available  to 
provide  output.  Most  UNIX  systems  provide  it. 

11.14.1  Temperature  Conversion 

/*.  ftoc  -  fahrenheit  to  Celsius  conversion,  -40f  to  lOOf  */ 

main( ) 

{ 

float  f,c 
f  =  -40.0; 

while  (  f  <=  100.0  )  { 
c  =  5.0*(f  -  32.0)/9.0; 
printf ("*4. Of  *6. 1f\n" ,f ,e) ; 
f  =  f  +  5.0; 

} 


.N  i 


.> 


c 


Examples 


£ 


11.1*1.2  Nearest  Points  In  this  example,  please  note  that  you  must  provide  the 
function  'gtarray'  (to  read  the  values  into  the  x  array),  and  that  the  'sqrt' 
function  must  be  made  available  to  your  program  by  accessing  the  appropriate 
system  library.  C  doesn't  have  an  exponentiation  operator,  either,  so  you  must 
either  multiply  the  terms  yourself  (as  I've  done  here)  or  provide  a  function  to 
do  it. 

/»  near  -  find  2  of  10  points  which  are  closest  neighbors  */ 
main  () 

{ 

float  xT 1 ol r 31 ;  /*  array  of  points  */ 

float  oldd,newd;  /*  distances  between  points  */ 

int  i,j,k,m,n;  /*  counters  */ 

gtarray(x, 10,3) ;  /*  function  to  fill  array  'x'; 

vou  must  provide  this  per 
your  operating  system  reqts.  */ 

/*  first  distance:  from  1  to  2  */ 

oldd  =  0.0; 

for^k  =  1;  k  <=  3;  ++k  ) 

oldd  =  oldd  +  C xT nfkl-xr2l[kl)*(x[ IITkl-xf ?lfkl ) ; 
oldd  =  3qrt(oldd); 

for(i  =  1;  i  <=  0;  ++i  )  ( 

ford  =  1  +  1;  .1  <=  10;  ++1  )  f 

newd  =  0.0; 

for ( k  =  1;  k  <=  3;  ++l<  ) 

newd  =  newd  +  ( x[i 1 [k l-*r 1 If k] ) *( xT il [k l-x(  |  IT k 1 ) ; 
newd  =  sqrt (newd); 

L f (  newd  <=  oldd  )  {  /*  found  a  new  nearest  neighbor  */ 

m  =  i ; 
n  =  .1; 

oldd  =  newd; 

1 

} 

} 

printf ("  closest  points  are  at  \n"); 

printf  ("id  id  1td\n",x[mir  n,x[ml[2l,x[mir3‘>); 

printf  ("id  id  id\n" ,x[nl[ 11 ,x[nl[ 2 ] , x[nlT 31 ) ; 

) 
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11.14.3  Count  the  'A 's 

/*.  counta  -  count  the  'A 'a  in  input.  •/ 
main( ) 

{ 

int  c,na; 
na  =  0; 

while(  (csgetcharO)  !  =  EOF  ) 
if(  c  ==  'A'  ) 

++na; 

printf("#d\n",na) 

} 

where  'getehar'  is  a  function  which  gets  the  next  character  from  <input>  and 
puts  it  into  the  variable  'e.'  If  <input>  is  at  end  of  file,  then  'getehar* 
assigns  a  non-character  value  to  e  (a  system  dependent  constant  of  0  or  -1).  I 
avoid  this  issue  by  using  the  symbolic  EOF. 


v.vcv-  >>>'.- 
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12.  FORTRAN77 


Blocks  of  statements  in  F0RTRAN77  are  grouped  only  in  the  'if  ...  then  else' 
statement,  since  blocks  are  a  new  concept  in  FORTRAN.  It  is  advisable  to 
'create'  blocks  of  statements  by  using  the  FORTRAN77  do-nothing  CONTINUE  state¬ 
ment. 

F0RTRAN77  requires  that  the  text  of  a  program  be  contained  in  columns  7  through 
72.  That  is,  you  may  begin  your  code  in  any  column  after  column  6.  It  is  a 
good  idea  to  indent  your  code  so  that  the  logical  structure  of  the  program  is 
reflected  by  the  layout  of  the  text  on  a  page. 

Comments  are  provided  in  a  FORTRAN77  program  by  putting  a  'C'  or  **'  in 
column  1 . 

•Since  all  The  Constructs  are  not  available  directly  in  FORTRAN77,  you  may  find 
it  easier  to  use  ratfor  than  this  language. 

12.1  Variable  Names 

Names  of  variables  are  limited  to  a  maximum  of  six  characters.  Every  variable 
should  be  declared  to  be  of  a  specific  type  and  its  purpose  defined.  This  may 
be  done  easily  through  the  type  definitions  REAL,  INTEGER,  LOGICAL,  CHARACTER, 
COMPLEX,  and  DOUBLEPRECIS ION.  It  is  inadvisable  to  allow  FORTRAN  to  type  your 
variables. 

The  type  definitions  may  al30  be  used  to  declare  the  dimensions  of  arrays,  so 
the  DIMENSION  statement  is  not  needed  in  FORTRAN  programs  (aee,  for  example,  the 
program  in  Section  12.14.2). 


Boolean  Expression 
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12.2  Boolean  Expression 
X  and  Y  are  variables  of  the  same  type. 
X  is  greater  than  Ys 
X  .GT.  Y 

X  is  greater  than  or  equal  to  Y: 

X  .GE.  Y 
X  is  equal  to  Ys 
X  .EQ.  Y 

X  is  not  equal  to  Ys 
X  .NE.  Y 

X  is  less  than  Y: 

X  .LT.  Y 

X  is  less  than  or  equal  to  Y: 

X  .LE.  Y 
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F0RTRAN77  Assignment  Statement 

12.3  Assignment  Statement 
X  =  expression 

FORTRAN77  statements  end  at  the  end  of  a  line.  Long  statements  may  be  continued 
for  up  to  19  additional  lines  by  providing  a  non-blank,  non-zero  character  in 
column  6  of  the  continuation  lines  (columns  1  through  5  must  be  blank).  Note 
however,  that  long  F0RTRAN77  statements  are  like  long  sentences:  they're  hard  to 
read,  hard  to  understand,  and  hard  to  correct!  Strive  for  short  statements, 
even  if  it  means  inventing  new  variables  to  contain  the  results  of  intermediate 
calculations. 

12.4  Conditional 

if(  B  )  then 

51 

52 
endif 

12.5  Alternative 

if  (  B  )  then 

51 

52 

else 

53 

54 
endif 

12.6  Multiple  Choice 

if  (  B1  )  then 

51 

else  if  (  B2  )  then 

52 

else  if  (  B3  )  then 

53 
•  •  • 

else 

Sn 

endif 
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12.7  While 

This  construct  must  be  simulated.  Because  a. reader  may  not  immediately  recog¬ 
nize  the  construct,  it  should  be  annotated  as  a  While.  Note  that  F0RTRAN77 

forces  you  to  express  the  Constructs  in  terms  of  transportation  -  "GOTO  label" 
is  a  rather  indirect  way  to  express  a  logical  concept. 

C  while  B 

23000  if  (  B  )  then 

51 

52 

goto  23000 
endif 

C  end  while 

12.8  Repea t-Ontil 

This  construct  must  be  simulated.  Because  a  reader  may  not  immediately  recog¬ 
nize  the  construct,  it  should  be  annotated  as  a  Repeat-Ontil . 

C  repeat  until  B 

23000  continue 

51 

52 

if  (  B  )  goto  23000 

12.9  Determinate  Iteration 

do  23000  I  =  J,K,M 

51 

52 

23000  continue 

J  is  the  initial  value  of  I,  K  is  the  terminal  value  of  I,  and  M  is  the  non-zero 
increment  of  I.  Note  that  if  the  initial  condition  satisfies  the  terminal  con¬ 
dition,  the  loop  will  not  be  executed. 
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F0RTRAN77  Altered  Loops 

12.10  Altered  Loops 
12.10.1  Loop  Exit 

Although  any  loop  (While,  Repeat-Until,  Determinate)  may  be  'exit'ed,  we  illus¬ 
trate  with  a  While. 


! 

23000  if  (  B1 
SI 

)  then 

\. 

if  ( 

B2  )  goto  23001 

>. 

S2 

goto 

23000 

. 

u’ 

endif 

23001  continue 
12.10.2  Loop  Redo 

Although  any  loop  (While,  Repeat-tint il ,  Determinate)  may  be  ' re-do 'ed,  we  illus¬ 
trate  with  a  While. 

23000  if  (  B1  )  then 

51 

if  (  B2  )  goto  23000 

52 

goto  23000 
endtf 
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Conditional:  if  (B)  then 

S 

endif 

Alternative:  if  (B)  then 

51 
else 

52 
endif 

Multiple  Choice:  if  (  B1  )  then 

51 

else  if  (  B2  )  then 

52 

else  if  (  B3  )  then 

53 

e  •  e 

else 

Sn 

endif 

While:  23000  if  (  B  )  then 

51 

52 

goto  23000 
endif 

Repeat-Until:  23000  continue 

51 

52 

if  (  B  )  goto  23000 

Determinate:  do  23000  I  =  J,K,M 

51 

52 

23000  continue 
Loop  Exit:  goto  label 

Loop  Redo:  goto  label 


FORTRAN77 


Program  Header 


12.12  Program  Header 

The  program  header  on  CDC  machines  has  the  form 

program  narae(  filel,  file2,  ...  ) 
declaration  of  variable  types 

12.13  Running  a  F0RTRAN77  Program 

Assume  you  have  composed  your  program  with  a  text  editor  on  a  CDC  computer  and 
the  code  is  in  <p>.  Assume  further  that  you  have  already  executed  the  commands 

ATTACH , ROTH , CCLLIB , ID=CSPR . 

LIBRARY, ROTH. 

(you  must  execute  these  commands  only  once).  Then  the  program  in  <p>  may  be  run 
with  the  command 

vc,p. 

abs. 


(see  [Roth83l) 
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12.14.1  Temperature  Conversion 


PROGRAM  FTOC 


C.  FTOC  -  FAHRENHEIT  TO  CELSIUS  CONVERSION,  -40F  TO  100F 
C 


REAL  F,C 


F  =  FAHRENHEIT,  C  s  CELSIUS 


F  =  -40.0 

IF(F  .LE.  100.0)  THEN 
C  =  5.0»(F  -  32.0)/9.0 
PRINT  *,F,C 
F  =  F  +  5.0 
GOTO  1 
ENDIF 


WHILE(  F  <=  100  ) 


END  WHILE 
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12.14.2  Nearest  Points  Assume  that  the  three  coordinates  of  each 
typed  on  a  single  line  of  <input>. 

PROGRAM  NEAR(  INPUT,  OUTPUT  ) 

C 

C.  NEAR  -  FIND  2  OF  10  POINTS  WHICH  ARE  CLOSEST  NEIGHBORS 
C 

C  X( 10,3)  IS  THE  ARRAY  OF  POINTS 
C  OLDD, NEWD  ARE  DISTANCES  BETWEEN  POINTS 
C  I,J,K,M,N  ARE  COUNTERS 
C 

REAL  X( 10,3) 

REAL  OLDD, NEWD 
INTEGER  I, J,K,M,N 
C 

READ  *,((X(I,J),J=1,3),I=1,10) 

C  FIRST  DISTANCE:  1  TO  2 

OLDD  =  SQRT(  (X( 1 , 1 )-X(2, 1 ) )»»2 
$  +  (X( 1 , 2)-X( 2, 2) )**2 

$  +  (X( 1,3)-X(2,3))**2  ) 

C 

DO  2  1=1,9 

IP1  =1+1 
DO  1  J  =  IP1,  10 

NEWD  =  SQRT(  (X( I , 1 )-X( J, 1 ) ) **2 

$  +  (X(I,2)-X(J,2))**2 

$  +  (X(I,3)-X(J,3))**2  ) 

c  NEW  MINIMUM  FOUND? 

IF(  NEWD  .LE.  OLDD  )  THEN 
M  =  I 
N  =  J 

OLDD  =  NEWD 
ENDIF 

1  CONTINUE 

2  CONTINUE 

PRINT  CLOSEST  POINTS  ARE  AT  ' 

PRINT  *,  fx(M,I),I=1,3),  (X(N,I),T=1 ,3) 

END 
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12. 1H. 3  Count  the  'A's  Approach:  assume  that  each  line  has  at  most  132  charac¬ 
ters,  and  read  an  entire  line  at  a  time.  Also  assume  that  character  variables 
can  be  represented  as  integers  (this  works  on  CDC  machines).  If  the  line  is 
shorter  than  132  characters,  FORTRAN  will  act  as  if  the  nonexistent  characters 
were  blanks.  If  the  line  is  longer  than  132  characters,  the  excess  will  be 
(silently)  truncated  (perhaps  not  the  best  solution). 

PROGRAM  COUNT* 


u 


C.  COUNT  THE  NUMBER  OF  'A'S  IN  <INPUT>. 
C 

INTEGER  COUNT, I 
CHARACTER  LINE  »1 32 


COUNT  =  0 

OPEN (UN IT=5, FILE= 'INPUT ' ) 

WHILE  NOT  END  OF  FILE 

READ (5, 2, END= 3)  LINE 
FORMAT(A132) 

DO  4  I  s  1,  132 

IF(LINE(I:I)  .EQ.  'A')  THEN 
COUNT  s  COUNT  +  1 
ENDIF 
CONTINUE 
GOTO  1 

END  WHILE 

CONTINUE 
PRINT  »,  COUNT 
END 


13.  FORTRAN66 


If  there  is  any  way  to  avoid  writing  programs  in  this  language,  you  should  take 
it.  One  way  which  comes  readily  to  mind  is  the  alternative  'ratfor,'  which  pro¬ 
duces  FORTRAN66  but  allows  you  to  express  your  program  in  terms  more  nearly  like 
English.  Note  that  the  straight-forward  application  of  FORTRAN66  costs  you  time 
and  effort  because  you  must  convolute  your  natural  thought  processes  when  you 
write  the  program,  and  then  again  when  you  go  looking  for  errors. 

The  concept  of  'blocks  of  statements'  does  not  exist  in  FORTRAN66.  It  is  advis¬ 
able  to  'create'  blocks  of  statements  by  using  the  FORTRAN66  do-nothing  CONTINUE 
statement. 

FORTRAN66  requires  that  the  text  of  a  program  be  contained  in  columns  7  through 
72.  That  is,  you  may  begin  your  code  in  any  column  after  column  6.  It  is  a 
good  idea  to  Indent  your  code  so  that  the  logical  structure  of  the  program  is 
reflected  by  the  layout  of  text  on  a  oage. 

Comments  are  provided  in  a  FORTRAN66  program  by  putting  a  'C'  in  column  1. 

13.1  Variable  Names 

Names  of  variables  are  limited  to  a  maximum  of  six  characters.  Every  variable 
should  be  declared  to  be  of  a  specific  type  and  its  purpose  defined.  This  may 
be  done  easily  through  the  type  definitions  REAL,  INTEGER,  LOGICAL,  COMPLEX, 
DOUBLEPRECISION.  It  is  inadvisable  to  let  FORTRAN  tvpe  your  variables  for  vou. 

The  type  definitions  may  also  be  used  to  declare  the  dimensions  of  arrays,  so 
the  DIMENSION  statement  is  not  needed  in  FORTRAN  programs  (see,  for  example,  the 
program  in  section  13.14.2). 
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13.2  Boolean  Expression 
X  and  Y  are  variables  of  the  same  type. 

X  is  greater  than  Y: 

X  .GT.  Y 

X  is  greater  than  or  equal  to  Y: 

X  .QE.  Y 
X  is  equal  to  Y: 

X  .EQ.  Y 

X  is  not  equal  to  Y: 

X  .NE.  Y 

X  is  less  than  Ys 
X  .LT.  Y 

X  is  less  than  or  equal  to  Y: 

X  .LE.  Y 

13.3  Assignment  Statement 
X  =  expression 

FORTRAN66  statements  end  at  the  end  of  a  line.  Long  statements  may  be  continued 
for  up  to  19  additional  lines  by  providing  a  non-blank,  non-zero  character  In 
column  6  of  the  continuation  lines  {colisnns  1  through  5  must  be  blank).  Note 
however,  that  long  FORTRAN66  statements  are  like  long  sentences:  they're  hard  to 
read,  hard  to  understand,  and  hard  to  correct!  Strive  for  short  statements, 
even  if  it  means  inventing  new  variables  to  contain  the  results  of  intermediate 
calculations. 
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13.4  Conditional 

In  general,  FORTRAN66  requires  a  'reversal  of  thought'  when  expressing  condi¬ 
tional  tests  because  the  only  control  available  to  group  statements  is  the  GOTO. 
Hence,  to  control  several  statements,  the  expression  'if  B  then'  is  expressed 
'if  NOT  B  goto'.  The  language  also  forces  you  to  express  logic  in  terms  of 
transportation  -  "GOTO  label"  is  a  rather  indirect  way  to  express  a  logical  con¬ 
cept.  Compare  the  following  with  the  definition  of  the  Conditional  Construct  in 
Section  4.4  and  in  some  of  the  other  languages:  sections  9.4,  10.4,  and  11.4. 

if ( .not. (B) )goto  23000 

51 

52 

23000  continue 

13.5  Alternative 

if ( .not. (B) )goto  23000 
C  B  IS  TRUE 

51 

52 

goto  23001 

23000  continue 

C  B  IS  FALSE 

53 

54 

23001  continue 

Commentary  in  FORTRAN66  are  lines  with  the  symbol  'C'  in  column  1.  Because  the 
control  structures  are  more  difficult  to  understand  in  this  language,  they 
should  be  more  thoroughly  commented  than  in  other  languages. 

13*6  Multiple  Choice 

This  form  of  multiple  choice  makes  you  assign  a  statement  label  and  a  condition 
number  to  any  of  the  choices  you  wish  to  make. 

goto  (  1,2,3,. *.,n  ),  I 

1  SI 
goto  m 

2  S2 
goto  m 

3  S3 
goto  m 

•  •  • 

n  Sn 

m  continue 
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13.7  While 

This  construct  must  be  simulated.  Beoause  someone  who  reads  your  program  may 
not  immediately  recognize  the  construct,  it  should  be  annotated  as  a  While. 

C  WHILE  B  DO 

23000  if( .not. (B) )goto  23001 

51 

52 

goto  23000 

23001  continue 

13.8  Repeat-Until 

This  construct  must  be  simulated.  Because  a  reader  may  not  immediately  recog¬ 
nize  the  construct,  it  should  be  annotated  as  a  Repeat-Until. 

C 

23000 


23001 


REPEAT  UNTIL  (B) 

continue 

51 

52 

if(  .not.(B)  )goto  23000 
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13.9  Determinate  Iteration 


13.9.1  Increasing  Index 


do  23000  I  =  J,K,M 

51 

52 

23000  continue 

J  is  the  initial  value  of  I,  K  is  the  terminal  value  of  I,  and  M  is  the  positive 
increment  of  I.  Mote  that  this  loop  will  be  executed  at  least  once  (like  the 
Repeat-Until) , 

13.9.2  Decreasing  Index  A  severe  lack  in  FORTRAN66,  this  may  be  simulated  with 
the  While  loop.  Old  FORTRANers  may  prefer  to  simulate  this  with  a  Repeat-Until 
so  that  the  loop  executes  at  least  once  and  is  an  exact  complement  of  the  Do. 


LOOP  FROM  J  D0WNT0  It  BY  M 


23000  if(I.lt.K)goto  23002 

51 

52 

I=I-M 

goto  23000 
23002  continue 


\>V jr'*  V*» 


Altered  Loops 
13.10  Altered  Loopa 


FORTRAN66 


13.10.1  Loop  Exit 

Although  any  loop  (While,  Repeat-Until,  Determinate)  may  be  'exit'ed,*  we  illus¬ 
trate  with  a  While. 

23000  if ( .not. (B1 ) )goto  23001 

51 

if (  B2  )  goto  23001 

52 

goto  23000 

23001  continue 

13.10.2  Loop  Redo 

Although  any  loop  (While,  Repeat-Until,  Determinate)  may  be  're-do'ed,  we  illus¬ 
trate  with  a  While. 

23000  if ( .not. (B1 ) )goto  23001 

51 

if(  B2  )  goto  23000 

52 

goto  23000 

23001  continue 


•  You  can  GOTO  almost  anywhere  in  a  FORTRAN66  program! 
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13.11  FORTRAN66  Summary 


Conditional:  if(.not.(B))goto  23000 

51 

52 

23000  continue 

Alternative:  if( .not.(B))goto  23000 

51 

52 

goto  23001 

23000  continue 

53 
SM 

23001  continue 

Multiple  Choice:  goto  (  1,2,3»...,n  ),  I 

1  SI 
goto  m 

2  S2 
goto  m 

3  S3 
goto  m 

•  •  • 

n  Sn 

m  continue 

While:  23000  if(.not.(B))goto  23001 

51 

52 

goto  23000 
23001  continue 

Repeat-Until:  23000  continue 

51 

52 

if  (  .not.(B)  )  goto  23000 

Determinate:  do  23000  I  =  J,K,M 

51 

52 

23000  continue 
Loop  Exit:  goto  label 


Loop  Redo 


goto  label 
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13.12  Program  Header 

The  FORTRAN 6 6  program  header  on  CDC  machines. has  the  form 

program  name(  fllel,  file2,  ...  ) 
declaration  of  variable  types 

13.13  Running  a  FORTRAN66  Program 

Assume  you  have  composed  your  program  with  a  text  editor  on  a  CDC  computer  and 
the  code  is  in  <p>.  Assume  further  that  you  have  already  executed  the  commands 

ATTACH , ROTH , CCLLIB , IDsCSPR . 

LIBRARY, ROTH. 

(you  must  execute  these  commands  only  once).  Then  the  program  in  <p>  may  be  run 
with  the  command 


fc,p. 

abs. 

(see  [Roth83l). 
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13.14  Examples 

13-14.1  Temperature  Conversion 
PROGRAM  FTOC( OUTPUT) 

.  FTOC  -  FAHRENHEIT  TO  CELSIUS  CONVERSION,  -40F  TO  100F 

F  =  FAHRENHEIT,  C  =  CELSIUS 

REAL  F,C 
F  r  -40.0 

WHILE  F  <=  100 

IF(F  .GT.  100.0)  GOTO  2 
C  =  5.0*(F  -  32.0)/9 
PRINT  *,F,C 
F  r  F  +  5.0 
GOTO  1 

END  WHILE 

CONTINUE 
END 
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13.14.2  Nearest  Points  Assume  that  the  three  coordinates  of  each  point  are 
typed  on  a  single  line  of  <input>. 

PROGRAM  NEAR(  INPUT,  OUTPUT  ) 

.  NEAR  -  FIND  2  OF  10  POINTS  WHICH  ARE  CLOSEST  NEIGHBORS 

X( 10,3)  IS  THE  ARRAY  OF  POINTS 
0LDD,NEWD  ARE  DISTANCES  BETWEEN  POINTS 
I, J,K,M,N  ARE  COUNTERS 

REAL  X( 10,3) 

REAL  OLDD, NEWD 
INTEGER  I, J,K,M,N 

READ  *,((X(I,J),J=1,3),I=1,10) 

C  FIRST  DISTANCE  FROM  1  TO  2 

OLDD  =  SQRT(  (X( 1 , 1 )-X(2, 1 ) )»*2 

$  +  (X(1,2)-X(2,2) )**2 

$  +  (X(1,3)-XC2,3))**2  ) 

C 

DO  2  1=1,9 

IP1  =1+1 
DO  1  J  =  IP1,  10 

NEWD  =  SQRT(  (X(I,1)-X(J, 1) )**2 

$  +  (X(I,2)-X( J,2) )**2 

$  +  (X(I,3)-X( J,3))**2  ) 

C  NEW  MINIMUM  FOUND? 

IF( .NOT. (  NEWD  .LE.  OLDD  ))  GOTO  3 
M  =  I 
N  =  J 

OLDD  =  NEWD 
3  CONTINUE 

1  CONTINUE 

2  CONTINUE 

PRINT  CLOSEST  POINTS  ARE  AT  ' 

PRINT  •,  (X(M,I),I=1 ,3),  (X(N,I) ,1=1 ,3) 
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13.14.3  Count  the  A'sj  Approach:  assume  that  each  line  has  at  most  132  charac¬ 
ters,  and  read  an  entire  line  at  a  time.  Aa.so  assume  that  character  variables 
can  be  represented  as  Integers  (this  works  on  CDC  machines).  If  the  line  is 
shorter  than  132  characters,  FORTRAN  will  act  as  if  the  nonexistent  characters 
were  blanks.  If  the  line  is  longer  than  132  characters,  the  excess  will  be 
(silently)  truncated  (perhaps  not  the  best  solution). 

PROGRAM  COUNTA( INPUT, OUTPUT, TAPE5=INPUT) 

C 

C.  COUNT  THE  NUMBER  OF  'A'S  IN  <INPUT>. 

C 

INTEGER  COUNT,  I,  LINE(132) 

INTEGER  EOF 
C 

COUNT  =  0 

READ(5, 1 )  (LINE(I),I=1,132) 

1  F0RMAT(1??!\1) 

C  WHILE  NOT  END  OF  FILE 

2  IF(  EOF( 5)  .NE.  0)  GOTO  4 

DO  3  I  =  1,  132 

IF(LINE(I)  .EQ.  'A')  COUNT  =  COUNT  +  1 

3  CONTINUE 
READ(5, 1 )  (LINE(I) ,1=1 ,132) 

GOTO  2 

C  END  WHILE 

4  CONTINUE 
PRINT  »,  COUNT 
END 

where  we  note  that  the  integer  function  EOF  is  a  CDCism. 
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14.  BASIC 

There  are  many  versions  of  BASIC.  The  following  remarks  therefore  apply  to 
BASIC  in  a  general  sense;  for  specifies,  you  will  need  the  RASIC  manual  for  your 
system. 

The  concept  of  'blocks  of  statements'  doesn't  exist  in  BASIC.  Most  BASICS 
prohibit  the  indentation  of  code  to  reflect  the  logical  structure  of  a  program; 
thus  BASIC  programs  are  usually  difficult  to  read  and  write. 

Comments  are  provided  in  a  BASIC  program  by  putting  the  string  'REM'  (for 
REMark)  as  the  first  entry  following  the  line  number.  Because  the  control 
structures  are  more  difficult  to  understand  in  this  language,  they  should  be 
more  thoroughly  commented  than  in  other  languages.* 

Since  few  of  the  Constructs  are  available  directly  in  BASIC,  you  may  find  it 
easier  to  use  any  other  language. 

14.1  Variable  Names 

The  names  of  BASIC  variables  may  be  one  or  two  characters  long.  The  first  char¬ 
acter  must  be  a  letter,  and  the  second  may  be  a  letter  or  a  number.  Most  BASICS 
assume  all  variables  are  of  type  REAL  (i.e,  they  have  fractional  parts).  Char¬ 
acter  type  variables  usually  include  the  '$'  as  the  last  character  in  the  name 
'e.g.,  N$) . 


*  Unfortunately,  BASIC  is  often  run  on  machines  whose  memory  is  too  small  to 
hold  both  code  and  comments l 
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14.2  Boolean  Expression 
X  and  Y  are  variables  of  the  same  type. 
X  is  greater  than  Y: 

X  >  Y 

X  is  greater  than  or  equal  to  Y; 

X  >=  Y 

X  is  equal  to  Y: 

X  =  Y 

X  is  not  equal  to  Y: 

X  <>  Y 

X  is  less  than  Ys 
X  <  Y 

X  is  less  than  or  equal  to  Y: 


Assignment  Statement 

14.3  Assignment  Statement 
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LET  X  =  expression 

BASIC  statements  end  at  the  end-of-line.  Continuation  of  statements  is  often 
impossible,  so  you  need  to  simplify  expressions  by  inventing  new  variables  to 
hold  the  results  of  intermediate  oaleulations. 

14.4  Conditional 

In  general,  BASIC  requires  a  'reversal  of  thought'  when  expressing  conditional 
tests.  This  is  because  the  only  command  available  to  group  statements  is  the 
'{goto}  linenumber'.  In  many  implementations,  the  'goto'  is  not  supplied,  but 
understood.  Hence,  to  control  several  statements,  the  expression  'if  B  then'  is 
expressed  'if  NOT  B  linenumber'. 

120  if  not  B  then  150 
130  SI 
140  S2 

150  rem  ...continue 

14.5  Alternative 


10  if  not  B  then  50 
15  rem  B  IS  TRUE 

20  SI 
30  S2 

40  goto  80 

50  rem  B  IS  FALSE 

60  S3 
70  S4 

80  rem  continue 
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14.6  Multiple  Choice 

The  BASIC  multiple  choice  forces  you  to  assign  a  program  line  msaber  (in  the 
proper  numerical  sequence)  which  corresponds  to  the  condition  number  of  the 
choice  you  wish  to  make.  This  may  be  difficult  to  do  when  you  are  designing  a 
program  because  you  can  'run  out'  of  line  numbers.  Allow  'enough'  numbers. 

10  on  I  goto  (  20,40,60,90  ) 

20  SI 

30  goto  110 
40  S2 

50  goto  110 
60  S3 

70  goto  110 
90  S4 

100  goto  110 

110  rem  ...continue 

14.7  While 

This  construct  must  be  simulated. 

5  rem  WHILE  B  DO 

10  if  not  B  then  50 
20  SI 
30  S2 

40  goto  10 
50  rem  continue 

14.8  Repeat-Until 

This  construct  must  be  simulated. 

10  rem  REPEAT  UNTIL  B 

20  SI 
30  S2 

40  if  not  B  then  10 

14.9  Determinate  Iteration 

This  construct  is  handled  nicely  in  many  BASICS. 

10  for  I  =  J  to  K  step  M 
20  SI 
30  S2 
40  next  I 

J  is  the  initial  value  of  I,  K  is  the  terminal  value  of  I,  and  M  is  the  non-zero 
increment  of  I.  Note  that  if  the  initial  condition  satisfies  the  terminal  con¬ 
dition,  the  loop  will  not  be  executed. 
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Altered  Loops 
14.10  Altered  Loops 
14.10.1  Loop  Exits 


BASIC 


Although  any  loop  (While,  Repeat-Until ,  Determinate)  may  be  'exit'ed  (you  can 
GOTO  just  about  anywhere  in  a  BASIC  program) ,  we  illustrate  with  a  While. 

10  if  not  B1  then  60 
20  SI 

30  if  B2  then  60 
40  S2 

50  goto  10 
60  rem  continue 

14.10.2  Loop  Redo 

Although  any  loop  (While,  Repeat-Until,  Determinate)  may  be  're-do'ed,  we  illus¬ 
trate  with  a  While. 

10  if  not  B1  then  60 
20  SI 

30  if  B2  then  10 
40  S2 

50  goto  10 
60  rem  continue 


BASIC 

14.11  BASIC  Sm«"»»*y 


Conditional: 

10 

if  not  B 

then  30 

20 

S 

30 

•  •  • 

Alternative: 

10 

if  not  B 

then  40 

20 

SI 

30 

goto  50 

40 

S2 

50 

•  •  • 

Multiple  Choice: 

10 

on  I  goto  (  20,40,60, 

20 

SI 

30 

goto  m 

40 

S2 

50 

goto  m 

n 

•  •  • 

Sn 

m 

•  •  • 

While: 

10 

if  not  B 

then  40 

20 

S 

30 

goto  10 

40 

•  •  • 

Repeat-Until: 

10 

S 

20 

if  not  B 

then  10 

Determinate: 

10 

for  I  =  . 

J  to  K  step  M 

20 

S 

30 

next  I 

Loop  Exit: 

goto  linenumber 

Loop  Redo: 

goto  linenumber 

BASIC  Summary 


Program  Header 


BASIC 


14.12  Program  Header 

BASIC  programs  do  not  require  a  header  on  most  systems.  However,  it  is  best  to 
use  at  least  one  comment  line  to  identify  the  program. 

14.13  Running  a  BASIC  Program 

Please  refer  to  the  manual  for  the  system  you  have.  BASIC  programs  are  usually 
entered  from  a  keyboard  or  recalled  from  auxiliary  storage  (tape  or  floppy 
disk).  The  sequence  to  run  a  program  is  usually  something  like 


RUN  prograraname 


j 


Examples 
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14.14.2  Nearest  Points  Assume  that  the  three  coordinates  of  each  point  are 
contained  in  DATA  statements  (some  points. are  provided  below).  It  is  usually 
easier  to  use  the  editing  capabilities  of  the  BASIC  interpreter  by  supplying 
DATA  statements  than  it  is  to  correctly  type  in  30  numbers  'interactively.' 


10  REM  NEAR  -  FIND  2  OF  10  POINTS  WHICH  ARE  CLOSEST  NEIGHBORS 
20  REM  X( 10,3)  IS  THE  ARRAY  OF  POINTS 
30  REM  C,D  ARE  DISTANCES  BETWEEN  POINTS 
40  REM  I,J,K,M,N  ARE  COUNTERS 
50  DIM  X( 10,3) 

65  REM  ...  GET  POINTS  FROM  DATA  STATEMENTS,  LINES  280-370 
70  FOR  I  =  1  TO  10 
80  FOR  J  =  1  TO  3 
90  READ  X(I,J) 

100  NEXT  J 


110  NEXT  I 

111  REM  .. 

112  C  s  0 

113  FOR  K 


FIRST  DISTANCE  FROM  1  TO  2 


113  FOR  K  =  1  TO  3 

114  C  =  C  +  (X( 1 ,K)-X(2,K) )“2 

115  NEXT  K 

1 16  C  =  SQRT(C) 

118  REM  ...  BEGIN  LOOKING  AT  NEIGHBORS 
120  FOR  I  =  1  TO  9 


118  REM  . 
120  FOR  I 
130  IP1  = 
140  FOR  J 
150  D  s  0 
160  FOR  K 


I  +  1 

r  IP1  TO  10 


160  FOR  K  r  1  TO  3 

170  D  =  D  +  (X(I,K)-X(J,K))*2 

180  NEXT  K 

190  D  =  SQRT(D) 

195  REM  ...  FOUND  NEW  MINIMUM? 
200  IF  D  >  C  THEN  240 


210  M  =  I 
220  N  =  J 
230  C  =  D 
240  NEXT  J 
245  REM  .. 
250  NEXT  I 
255  REM  .. 


END  FOR  J  =  IP1  TO  10 


255  REM  ...  END  FOR  I  =  1  TO  9 

260  PRINT  "  CLOSEST  POINTS  ARE  AT 

270  PRINT  X(M,1),X(M,2),X(M,3),X(N,1),X(N,2),X(N,3) 

280  DATA  10,20,30 

290  DATA  21,31,41 

300  DATA  32,42,52 

310  DATA  10.3,11,27 

320  DATA  -27.006,23,4 

330  DATA  22,9.3,26 

340  DATA  1,2, 2.1 

350  DATA  3,31,20 

360  DATA  3,31.1,-24 

370  DATA  19,-21.07,3.001 

380  END 
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14.14.3  Count  the  'A's  The  BASIC  language  recognizes  the  character  data  type. 
However,  the  language  does  not  have  a  standard  way  to  read  a  file;  data  are  usu¬ 
ally  embedded  in  the  program  itself .  To  access  files  on  a  system  running  BASIC, 
one  must  somehow  execute  calls  to  the  operating  system  through  'extensions'  to 
the  language.  Since  these  calls  vary  widely  from  vendor  to  vendor,  I  omit  this 
example. 
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17.  GLOSSARY 


algorithm  a  description  of  the  steps  necessary  to  perform  a  calcula¬ 

tion.  This  is  not  to  be  confused  with  algorasm,  which  is 
what  a  programmer  has  the  first  time  his  program  runs. 

code  programs,  or  fragments  of  programs,  are  called  'code'  because 

they  represent  ideas  in  a  form  which  can  be  read  by  a 
machine. 

compiler  a  program  which  translates  some  'high-level'  language  (such 

as  Pascal,  FORTRAN,  Ada)  into  the  equivalent  instructions’ 
that  a  computer  can  understand. 

computer  an  extremely  gullible  machine  which  often  does  exactly  what 

you  tell  it  to  do,  rather  than  what  you  really  want  it  to  do. 

down  the  state  in  which  a  computer  is  immune  to  user  input. 

English  an  obtuse  language  correctly  spoken  only  by  Edwin  Newman  and 

William  F  Buckley,  Jr. 

error  something  which  causes  a  program  to  not  compile  (sometimes 

called  a  typographical  error),  or  which  causes  a  program  to 
be  'not  running.'  Note  that  a  program  which  never  stops  is 
considered  to  be  a  'not  running'  program. 

execution  Jargon  for  'run.' 

function  (a.)  what  a  program  is  supposed  to  do.  (b.)  a  program  module 

(e.g.,  the  TANGENT  function). 

go  verb  describing  the  act  of  attempting  to  get  a  program  to 

run. 

implementation  how  a  program  is  put  on  a  particular  computer  system. 

interactive  In  an  'interactive'  computing  environment,  man  and  machine 

have  a  dialog:  man  types  command,  mach*-'  a  he  work,  man 
types  command,  etc.  As  opposed  to  'ba*  '  x  .\eh  man  per¬ 
forates  regularly  shaped  pieces,  of  pastebv.  .  d,  man  totes 
pasteboards  to  machine,  machine  looks  through  holes,  machine 
prints  what  it  mis-read,  man  totes  pasteboard  and  paper  back 
to  perforator  for  another  cycle. 

library  a  special  kind  of  file  in  which  one  may  store  programs,  sub¬ 

routines,  functions,  and  the  like.  Libraries  make  it  easy 
for  you  to  use  code  developed  by  others. 

line  (a.)  telephone  connection  to  a  computer,  (b.)  a  series  of 

characters  terminated  by  a  'line-feed'  character,  (c.)  the 
entire  contents  of  a  computer  punchcard. 
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load  a  CDC  verb.  To  'load'  a  program  is  to  put  a  machine  language 

program  into  memory  along  with  other  modules  from  other 
libraries.  On  other  systems,  this  is  known  as  'binding'  or 
'linking'  or  'link-editing.' 

loop  a  synonym  for  iteration. 

macro  a  fragment  of  code  which  is  'expanded'  by  a  program  called  a 

macro-processor.  Macros  are  used  to  extend  the  power  of  a 
language.  Simple  examples:  the  arithmetic  statement  function 
of  FORTRAN66  and  FORTRAN77,  and  the  PARAMETER  statement  of 
FORTRAN77. 

memory  is  the  place  where  the  oomputer  'stores'  data.  Because  read¬ 

ing  the  computer  memory  doesn't  destroy  what  is  read,  one  can 
perceive  that  the  computer  'remembers'  numbers.  Memory  i3 
often  called  'prime  store'  in  British  publications. 

operating  system  a  computer  program  which  'runs  the  machine'  and  allows  other 

programs  to  run. 

preprocessor  a  program  which  runs  before  the  'real'  work  starts.  For 

example,  a  text  editor  is  a  preprocessor  to  ratfor,  ratfor  is 
a  preprocessor  to  FORTRAN,  FORTRAN  is  a  preprocessor  to  the 
LOADER,  the  LOADER  is  a  preprocessor  to  the  SYSTEM,  the  SYS¬ 
TEM  is  a  preprocessor  (and  co-processor)  of  YOUR  PROGRAM. 

procedure  (a.)  a  Pascal  subroutine;  (b.)  a  CCL  program,  (c.)  a  way  of 

doing  things. 

program  a  set  of  definitions,  declarations,  data,  and  algorithms 

which  becomes  a  component  of  a  computer. 

recursion  the  expression  of  a  function  in  terms  of  itself.  E.g.,  the 

factorial  n!  may  be  expressed  recursively  as 

n!  =  n*(n-1)f  with  n  >  0  and  1!  a  1. 

Of  the  languages  discussed  here,  recursion  is  possible  only 
in  C  and  Pascal. 

register  a  computer  component  which  can  hold  a  data  item. 

robust  software  is  said  to  be  robust  when  it  can  withstand  the 

assault  of  any  data  without  'going  down.' 

run  Execution  of  a  program. 

running  Engineering  programs  are  'running'  when  they  produce  correct 

output.  Any  other  program  is  'not  running.'  See  also 
'error. ' 
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shell  a  program  which  provides  the  interactive  interface  between 

man  and  modern  computers. 

software  the  hard  part  of  computing  (as  opposed  to  hardware,  which  is 

the  easy  part). 

store  as  a  verb:  to  put  data  into  the  computer's  memory.  as  a 

noun:  the  place  where  data  are  kept.  See  also  'memory.' 

symbol  a  representation  of  a  single  thing;  a  symbol  may  be  composed 

of  more  than  one  character.  For  example,  the  Pascal  symbols 
used  to  group  statements  are  'begin'  and  'end.' 

system  the  computer  program  that  allows  other  computer  programs  to 

run.  Also  called  'executive'  and  'monitor.' 

word  a  grouping  of  information  on  a  machine  which  is  convenient  to 

the  manufacturer.  On  CDC  equipment,  a  word  is  ten  6-bit 
characters;  on  some  DEC  equipment,  a  word  is  four  8-bit 
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